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Flux  qubits:  controllable  coupling,  single-shot  readout  and  hot  electrons 

John  Clarke,1  T.  Hime,1,  S.  Linzen,1,  B.L.T.  Plourde,1,  P.A.  Reichardt,1,  T.L.  Robertson,1 
K.B.  Whaley,2,  F.K.  Wilhelm,3,  C.-E.  Wu,1,  and  J.  Zhang2 

Departments  of 1  Physics  and  2  Chemistry,  University  of  California,  Berkeley,  C A  94720,  USA 
3Ludwig-Maximillians-Universitdt,  80333  Miinchen,  Germany 

The  flux  state  of  a  three-junction  qubit  is  measured  with  a  dc  Superconducting  QUantum 
Interference  Device.  The  qubit  and  SQUID  are  fabricated  from  evaporated  aluminum  films 
that  are  patterned  using  electron-beam  lithography.  The  qubit  encloses  an  area  of  11280  pm  , 
so  that  on-chip  flux  lines  can  be  used  to  generate  the  necessary  half-flux  quantum  of  bias 
(Fig.  1);  nonetheless,  the  geometrical  inductance  of  the  qubit  is  small  compared  with  the 
kinetic  inductance  of  the  three  junctions.  The  flux  state  of  the  qubit  is  determined  by 
applying  current  pulses  to  the  SQUID  and  adjusting  the  amplitude  of  the  pulse  to  obtain  a 
50%  switching  probability.  We  discuss  two  aspects  of  the  experiment:  single-shot  readout  of 
the  qubit  with  a  SQUID  in  which  each  junction  is  shunted  with  a  frequency-dependent 
impedance,  and  the  role  of  hot  electrons  in  the  SQUID  in  reducing  the  relaxation  time  of  the 
qubit.  We  begin,  however,  with  a  proposed  scheme  for  varying  the  coupling  between  two 
flux  qubits  by  adjusting  the  bias  current  of  the  readout  SQUID  in  the  zero-voltage  state. 

Fig.  1:  SQUID  and  three-junction  flux  qubit. 
SQUID  bias  current  and  voltage  leads  at  top  and 
bottom  of  SQUID  loop.  Flux  bias  lines  located  at 
top  and  sides  of  SQUID  loop.  The  Josephson 
tunnel  junctions,  located  at  the  bottom  of  the 
SQUID  and  qubit  loops,  are  typically  200  x  200 
nm2  and  are  formed  with  double-angle 
evaporation. 


In  order  to  entangle  two  flux  qubits,  it  is  highly  desirable  to  be  able  to  control  their 
mutual  inductance  M.  In  particular,  manipulations  to  achieve  entanglement  are  much  more 
straightforward  if  M  can  be  reduced  precisely  to  zero  at  specified  times.  This  is  not  a  trivial 
problem  because  of  the  inherent  coupling  of  the  two  qubits  even  in  the  absence  of  external 
circuits.  A  scheme  that  enables  one  to  vary  the  coupling  and  to  reduce  it  to  zero  is  shown  in 
Fig.  2(a),  in  which  two  qubits  are  surrounded  by  a  SQUID.  The  technique  relies  on  the  fact 
that  the  inverse  dynamic  inductance  of  a  SQUID  in  the  zero- voltage  state  can  be  either 
positive  or  negative.  When  a  flux  is  applied  to  the  SQUID,  the  flux  generated  by  the  induced 
supercurrent  opposes  applied  flux  in  the  first  case,  but  supports  it  in  the  second.  Thus,  if 
there  is  a  change  in  the  current  in  (say)  qubit  1,  there  will  be  a  change  in  the  flux  in  qubit  2 
due  to  the  mutual  inductance  between  the  qubits,  and  a  second  contribution  that  is  mediated 
by  the  SQUID.  This  SQUID-mediated  contribution  can  be  made  positive  or  negative  simply 
by  changing  the  magnitude  of  the  bias  current  in  the  SQUID  at  an  appropriate,  fixed  value  of 
flux  [Fig.  2(b)];  the  SQUID  remains  in  the  zero-voltage  state.  Since  the  change  in  bias 
current  required  to  control  the  inductive  coupling  of  the  two  qubits  is  small,  less  than  the 
critical  current  of  the  SQUID,  it  can  be  pulsed  rapidly  with  existing  technology.  A  set  of 


device  parameters  is  proposed  that  is  realistically  achievable.  This  method  enables  one  both 
to  vary  the  mutual  inductance  between  two  qubits  and  to  read  their  flux  states  with  a  single 
SQUID,  simply  by  applying  current  pulses  of  appropriate  magnitude.  A  scheme  for  a 
Controlled  Not  (CNOT)  gate  is  suggested  that  involves  only  current  pulses  and  precisely 
chosen  microwave  pulses  (Fig.  3). 


Fig.  2:  (a)  SQUID-based  coupling  scheme.  The 
SQUID  has  inductance  L  and  critical  current  Ic. 
The  admittance  Y  represents  the  SQUID  bias 
circuitry,  (b)  Response  of  SQUID  circulating 
current  J  to  applied  flux  ®s  for  pL  =  2  L  Io  /  ®o  = 
0.04  and  Ib/Ic  (0.45  ®o) =  0,  0.4,  0.6,  0.85  (top  to 
bottom).  Lower  right  inset  shows  J(<t>s)  for  same 
values  of  lb  near  ®s  =  0.45  ®0-  Upper  left  inset 
shows  Ic  versus  ®s. 


Fig  .3  :  Pulse  sequence  for  implementing 
CNOT  gate.  Energy  scales  in  GHz.  The 
bias  current  is  pulsed  to  turn  on  the 
interaction  in  the  central  region.  Total 
single-qubit  energy  bias  £j(t)  =  E;°  +  £jm(t) 

+  8si(t),  where  e®  are  the  static  qubit  bias 
energies  and  £i,2m(t)  are  microwave  pulses 
which  produce  single-qubit  rotations  in  the 
decoupled  configuration;  crosstalk  shift  of 
Ej(t)  due  to  bias  current  pulse,  8£;(t),  along 
with  modulation  of  coupling  energy  K(t) 
due  to  microwaves  are  shown. 


The  use  of  a  dc  SQUID  with  unshunted  tunnel  junctions  has  the  advantage  of  creating 
low  dissipation  in  the  flux  qubit  but  the  disadvantage  of  a  broad  switching  distribution  due  to 
the  onset  of  macroscopic  quantum  tunneling  (MQT)  at  low  temperatures.  Adding  resistors 
across  the  junctions  narrows  the  distribution  but  substantially  increases  dissipation  in  the 
qubit.  We  demonstrate  that  a  resistor  R  in  series  with  a  capacitor  C  connected  across  each 
junction  enables  one  to  achieve  single-shot  readout  while  maintaining  low  dissipation  in  the 
qubit.  The  RC-time  constant  is  chosen  so  that  the  reactance  of  the  capacitor  is  small  at  the 
plasma  frequency  of  the  junctions,  ~100  GHz,  providing  substantial  damping  that  suppresses 
MQT.  On  the  other  hand,  at  the  qubit  frequency,  ~1  GHz,  the  reactance  of  the  capacitor  is 
much  greater  than  R,  and  the  dissipation  coupled  to  the  qubit  is  low.  Detailed  calculations  of 
the  width  of  the  switching  distribution  and  of  the  qubit  relaxation  and  dephasing  rates  due  to 
the  RC  shunts  are  presented. 

When  the  SQUID  switches  to  the  gap  voltage  2A/e  (A  is  the  energy  gap  of  the  A1  films) 
following  a  measurement  of  its  critical  current,  there  is  substantial  dissipation.  In 
spectroscopic  measurements  on  a  flux  qubit,  we  have  measured  the  relaxation  of  the 
microwave-induced  peaks  as  a  function  of  the  delay  time  t  between  successive  current  pulses 
in  the  readout  SQUID.  We  find  that  the  relaxation  rate  increases  dramatically  as  t  is  reduced. 
We  explain  this  behavior  in  terms  of  a  model  in  which  hot  electrons,  generated  when  the 
SQUID  switches  to  2A/e,  increase  the  subgap  conductance  of  the  tunnel  junctions.  This 
conductance  increases  the  dissipation  experienced  by  the  flux  qubit  for  a  remarkably  long 
time,  ~1  ms,  after  the  SQUID  has  been  restored  to  its  zero- voltage  state. 

This  work  was  supported  by  AFOSR,  ARO  and  NSF. 
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We  have  studied  quantum  coherent  dynamics  of  two  coupled  superconducting 
charge  qubits1.  Each  qubit  is  based  on  a  Cooper-pair  box,  a  small  superconducting 
island,  coupled  to  a  reservoir  through  a  small  Josephson  junction.  The  qubits  are 
electrostatically  coupled  to  each  other  by  a  miniature  on-chip  capacitor  formed  by 
overlapping  A1  layers.  The  sample  was  fabricated  by  an  angle  evaporation  of  A1  layers 
through  a  suspended  mask  formed  by  conventional  electron-beam  lithography. 

In  the  charge  regime  when  the  charging  energy  of  each  qubit  E\t2  exceeds  the 
Josephson  energyEji^,  the  Hamiltonian  of  the  system  can  be  reduced  to  the  four-state 
Hamiltonian: 

H  =  I  Enl>t2\nin2){n\n2 1 ~ ~ir  I(|O)(lM#|)0h)(”2|- 

«1, ”2=0,1  l  ”2=0,1  ^ 

%  Xh}(»iH0X1l+|i>(°|) 

where  Emm  is  the  electrostatic  energy  of  the  system  corresponding  to  the  charge 
configuration  (n\,n2),  and  the  four  charge  states  |00),  |10),  |0l}  and  |ll)  are  used  as  a 

basis.  Here  the  first  and  the  second  indices  refer  to  the  number  of  Cooper  pairs  in  the 
first  and  the  second  qubit,  respectively.  By  changing  dc  gate  voltages  we  can  control  the 
initial  state  of  the  system.  There  is  a  special  point  in  the  energy  diagram,  that  we  call  a 
double  degeneracy  point,  where  £oo  =  E\\  and  E\q  -  £oi-  By  preparing  the  system  in  the 
1 00}  initial  state  and  then  bringing  the  system  non-adiabatically  to  the  double 

degeneracy  point,  we  create  conditions  for  coherent  evolution  of  the  system.  Non- 
adiabatic  shift  of  the  system  is  done  by  a  sharp  control  pulse  applied  to  both  qubits 
simultaneously.  As  a  result  of  the  evolution,  the  final  state  of  the  system  is  a 
superposition  of  all  four  charge  states:  =  CjjOO)  c2jl0)  -f-  c3j0l)  -»-  c4jl  l) ,  where  the 

coefficients  c\  depend  on  time  and  hence  on  the  length  of  the  control  pulse.  The  problem 
can  be  solved  analytically  for  a  rectangular  pulse  shape.  To  readout  each  qubit,  we 
attach  a  properly  biased  probe  electrode  to  each  Cooper  pair  box  through  a  resistive 
tunnel  junction.  We  then  measure  pulse  induced  probe  currents  I\  and  /2  through  the 
probes.  Each  current  is  proportional  to  the  probability  _pi,2(l)  of  the  corresponding  qubit 

to  be  in  state  jl) ,  i.e.,  px  (1)  =  c\  (t)  +  c\  (t)  and  p2  (1)  =  c2  (t)  +  c\  (/) .  Thus,  by  measuring 
the  probe  currents  we  trace  pi>2(l). 


Eji  ^Hz) 


Fig.  1 .  Avoided  level  crossing  of  two  coupled  charge  qubits. 


The  observed  oscillations  of  I\  and  h  contained  two  frequency  components  in  the 
spectrum,  Q  +  e  and  Q  +  £,  that  depend  on  the  sample  parameters  and  result  in  beatings. 
By  tracing  oscillations  for  different  values  of  En,  we  observed  avoided  level  crossing  in 
the  energy  spectrum  as  shown  in  Fig.  1.  The  observed  beatings  and  avoided  level 
crossing  are  direct  evidence  for  the  interaction  between  the  qubits  and  in  a  remarkable 
quantitative  agreement  with  the  theoretical  expectations. 

Our  experimental  results  show  that  the  coherence  time  of  the  coupled  qubit 
oscillations  is  a  factor  of  four  shorter  compared  to  the  case  of  single  qubit  oscillations. 
This  fact  is  still  to  be  understood.  Nevertheless,  even  with  short  coherence  time  we  were 
able  to  perform  logic  operation. 

By  adding  one  more  pulse  gate  to  the  coupled  qubit  circuit  we  were  able  to  address 
each  qubit  individually.  This  allowed  us  to  demonstrate  conditional  gate  operation,  a 
prototype  of  the  quantum  C-NOT  logic  gate. 

Our  simulations1 2  show  that  in  the  course  of  oscillations  the  qubits  remain  entangled 
most  of  the  time. 


[1] .  Yu.  A.  Pashkin,  T.  Yamamoto,  O.  Astafiev,  Y.  Nakamura  et  al.  Nature  421,  823  (2003). 

[2] .  Yu.  A.  Pashkin,  T.  Tilma,  D.  V.  Averin,  O.  Astafiev,  T.  Yamamoto  et  al.  IJQI 1, 421  (2003). 
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A  variety  of  Josephson-junction-based  qubits  have  recently  been  implemented 
with  remarkable  progress  in  coherence  time  and  read-out  scheme.  These  developments, 
together  with  its  potential  scalability,  have  renewed  the  position  of  this  solid-state 
device  as  a  strong  candidate  as  a  building  block  for  the  quantum  computer.  On  the  other 
hand,  coupling  multiple  qubits  to  construct  a  logic  gate  is  another  important  step  toward 
the  realization  of  quantum  computer.  In  this  talk,  we  introduce  our  experiments  on 
conditional  gate  operations  using  coupled  superconducting  charge  qubits[l]. 

Samples  are  fabricated  by  standard  e-beam  lithography  and  shadow 
evaporation  technique.  A  scanning-electron  micrograph  of  the  device  is  shown  in  Fig.  1 . 
Blight  parts  are  made  of  Al.  Two  Cooper-pair  boxes  are  electrostatically  coupled  by  an 
on-chip  capacitor  which  is  shown  in  the  dashed  circle  in  the  figure.  The  right  qubit  has 
SQUID  geometry  and  we  use  this  qubit  as  the  control  qubit  and  the  left  one  as  the  target 
qubit.  Two  independent  pulse  gates  allow  us  to  address  each  qubit  individually.  The 
state  of  this  two-qubit  system  is  measured  by  currents  through  two  probes. 

The  logic  operation  we  tried  to  demonstrate  is  summarized  as  a  classical  truth 
table  in  Table  1 .  Target  qubit  is  flipped  only  when  the  control  qubit  is  in  the  |0>  state. 
This  operation  was  realized  by  utilizing  the  difference  of  the  degeneracy  condition 
between  two  pairs  of  the  charge  states,  namely,  |00>,  |01>  and  |10>,  |11>. 


Pulse  gate  2  Pulse  gate  1 

Figure  1 :  SEM  image  of  the  device. 
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Table  1 :  Truth  table  for  the  conditional 
gate. 


Figure  2(a)  shows  the  pulse  sequence,  which  we  utilized  in  the  experiment.  It 
consists  of  the  input  preparation  and  the  logic  operation.  The  first  pulse  applied  to  the 
Pulse  gate  1  creates  the  superposition  state  a|00>+p|10>.  This  state  should  be 
transformed  into  the  state  a|01>+(3|10>,  if  the  operation  pulse,  the  second  pulse  applied 
to  the  Pulse  gate  2,  is  described  by  Table  1 .  Figure  2(b)  shows  the  output  currents  of  the 
control  and  the  target  qubits  under  the  application  of  this  pulse  sequence  as  a  function  of 
the  Josephson  energy  of  the  control  qubit  En .  By  changing  En,  we  can  change  the 
coefficients  a  and  (3,  while  keeping  the  pulse  length  constant.  Observed  anti-correlation 
between  the  magnitudes  of  the  output  currents  is  consistent  with  the  above  expectation, 
indicating  that  our  operation  pulse  is  properly  working  as  a  conditional  gate. 
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Figure  2:  (a)  The  pulse  sequence  and  (b)  the  output  currents  of 
the  control  and  the  target  qubits. 


[1]  T.  Yamamoto,  Yu.  A.  Pashkin,  O.  Astafiev,  Y.  Nakamura  and  J.  S.  Tsai,  Nature  425, 
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Josephson  charge-phase  qubit  with  radio  frequency  readout: 
optimization  of  coupling  and  decoherence 
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The  two-junction-SQUID  configuration  of  the  Cooper-pair-box  qubit  enables  not  only 
the  control  of  Josephson  coupling,  but  also  a  readout  of  the  quantum  states  taking 
advantage  of  the  phase  degree  of  freedom.  This  readout  can  be  realized  by  measuring 
either  the  persistent  Josephson  current  circulating  in  the  superconducting  loop  [1]  or  the 
Josephson  inductance  of  the  qubit  junctions  with  small  island  in  between,  i.e.  the  single 
Cooper  pair  transistor,  closed  by  this  loop  [2],  The  Josephson  inductance  readout  is 
performed  by  inducing  small  radio  frequency  oscillations  of  frequency  ©  in  the  loop  and 
measuring  the  resonance  frequency  shift  in  a  weakly  coupled  tank  circuit.  Although 
such  measurement  requires  sufficiently  long  time,  presumably  »  Q/©,  where  Q  is  the 
quality  factor  of  the  tank,  the  rate  of  decoherence  caused  by  the  setup  is  significantly 
reduced.  This  is  achieved  by  appropriate  choice  of  the  operation  point  on  the  {q,  d>) 
plane,  where  q  is  the  gate  charge  and  <D  is  the  flux  applied  to  the  loop.  Specifically, 
operation  in  the  optimal  (magic)  points  [1,3]  makes  the  qubit  much  less  sensitive  to  the 
charge  and  phase  fluctuations  and  should  even  permit  the  single  shot  measurements. 
The  problems  of  asymmetry  in  the  junction  parameters,  quasiparticle  tunneling  and 
back  action  of  the  readout  setup  will  be  addressed  in  the  talk. 

[1]  D.  Vion,  A.  Aassime,  A.  Cottet,  P.  Joyez,  H.  Pothier,  C.  Urbina,  D.  Esteve  and  M.H.  Devoret, 

Science  296,  886  (2002). 

[2]  A.B.  Zorin,  Physica  C  368,  284  (2002). 

[3]  A.B.  Zorin,  arXive:cond-mat/03 12225  . 
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Fast  and  accurate  pumping  of  charge  in  Josephson  junction  circuits  poses  a  twofold 
problem.  On  one  hand  the  speed  is  typically  limited  by  the  magnitude  of  Josephson 
coupling  E}  due  to  Landau-Zener  crossing,  whereas  errors  due  to  supercurrent  leakage 
and  quantum  interference  effects  increase  with  increasing  Ej.  Moreover,  simultaneous 
suppression  of  the  latter  two  is  difficult  due  to  their  mutually  orthogonal  dependence 
on  Josephson  phase  difference.  To  circumvent  these  apparent  problems  we  propose 
and  test  an  alternative  device  concept,  the  Cooper  pair  sluice,  shown  schematically  in 
Fig.  1. 


Figure  1.  (a)  Basic  configuration  of  a  Cooper  pair  sluice,  (b)  Three  junction  SQUIDs 
for  better  suppression  of  Josephson  coupling,  (c)  Gate  and  flux  cycles  to  pump  one 
Cooper  pair  from  right  to  left  through  the  island. 

The  central  idea  of  the  sluice  is  to  vary  the  Josephson  couplings  in  time  in  order  to 
allow  fast  charge  transport  when  necessary  but  to  suppress  leakage  contributions  at 
other  times.  Therefore  the  sluice  consists  of  two  tunable  mesoscopic  junctions  (DC- 
SQUIDs),  with  just  one  island  in  between  [Fig.  1  (a)].  An  example  of  an  operating 
cycle  is  shown  in  Fig.  1  (c).  Initially  the  system  is  in  the  neutral  charge  state  with  ng 
(gate  charge)  and  Josephson  couplings  £jr  and  Eji  set  to  zero.  Thereafter  one  of  the 
Josephson  couplings  Ejr  is  tuned  adiabatically  to  maximum,  and  simultaneously  ng  is 
ramped  up  also.  In  the  middle  of  the  cycle  Eir  has  been  turned  off  again.  At  this 
moment  when  ESr  =  0  again,  exactly  one  pair  (or  in  general  an  integer  number  of  pairs 
depending  on  the  amplitude  of  nz)  has  been  brought  from  the  right  source  to  the 
central  island.  In  the  second  half  of  the  cycle  the  pair(s)  is  (are)  taken  out  through  the 
left  SQUID  in  a  similar  fashion. 

We  demonstrate  theoretically  that  it  is  possible  to  obtain  a  current  as  high  as  0.1  nA  at 
better  than  ppm  accuracy  using  a  sluice  (and  standard  aluminium  technology).  We 
discuss  the  influence  of  imperfect  suppression  of  the  Josephson  coupling  and  the  finite 
operating  frequency. 


Figure  2  shows  scanning  electron  micrograph  of  the  SQUIDs  and  input  coils  of  the 
sluice  in  (a)  and  a  blow-up  of  the  island  region  in  (b).  The  tickmark  in  (b)  is  1  pm 
long.  The  two-dimensional  flux-modulation  of  the  supercurrent  through  the  sluice  in 
(c)  indicates  better  than  90  %  suppression  of  Josephson  couplings  and  only  very  weak 
cross-coupling  between  the  SQUIDs  and  their  input  coils. 


Figure  2.  (a)  and  (b)  show  a  fabricated  sluice  structure.  Supercurrent  through  the 
sluice  plotted  against  the  two  applied  fluxes  (i.e.,  currents  in  the  input  coils). 

The  first  pumping  measurements  have  demonstrated  the  operation  of  the  device  but 
still  with  non-optimal  circuit  parameters. 
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Electrometers  based  on  Cooper  pair  tunneling  are,  in  principle,  dissipationless,  and  thus 
they  are  expected  to  exhibit  only  a  small  back  action  noise  [1],  We  have  studied  the 
dynamics  of  inductive  single  electron  transistor  (L-SET)  which  is  one  candidate  for  this 
kind  of  electrometry. 

The  L-SET  consists  of  a  superconducting  SET  (SSET)  connected  in  parallel  with  a  LC 
resonant  circuit  (See  Fig.  la).  The  SSET  is  described  by  the  Hamiltonian  is  ( q  -  qgf/2Cz 
-  2£jcos(cp/2)cos(0).  Here,  qg  is  the  gate  charge,  cp  is  the  phase  difference  across  the 
whole  device,  assumed  to  be  a  classical  variable  due  to  an  environment  having 
impedance  much  smaller  than  Rq  =  h/(2ef  ~  6.5  kO,  whereas  0,  conjugate  to  the  charge 
of  the  island  q,  experiences  quantum  effects.  The  lowest  of  the  resulting  bands  En(ty,qg) 
has  approximately  sinusoidal  energy  Eo(q>,qg)  —  Efqg)  cos(cp)  and  thus,  the  SSET  is 
effectively  a  gate-tunable  single  junction  [2],  From  our  point  of  view  the  most  important 
property  of  SSET  is  its  Josephson  inductance  Ljl((p,q^)  =  (2ti/Oo)2  c?E  ((p,qg)/d  (p2, 
where  <t>0  =  hl(2e)  is  the  flux  quantum.  With  a  total  shunting  capacitance  C,  a  SSET  thus 
forms  a  harmonic  oscillator  with  the  plasma  frequency /p  =  l/(2tt)  (Lj  C)"1/2  for  small  (p. 


FIGURE  1.  a)  Schematics  of  experimental  setup  used  in  our  measurements.  The  resonant  frequency  of  the 
LC  oscillator  is  613  MHz.  Coulomb  energy  of  the  SSET  Ec  =  0.92  K  and  the  ratio  E}  /  Ec  -  1.7.  b) 
Variation  of  the  phase  of  the  reflected  signal  from  the  L-SET  circuit  using  small,  harmonic  oscillation 
amplitudes  in  the  tank  circuit;  the  scale  of  variation  in  degrees  is  shown  on  the  right. 

Energy  of  the  system  of  the  SSET  plus  the  LC  oscillator  is  the  sum  of  potential  energy 
£n((p,<7g)  +  <I>2/(2L),  approximately  -  Efqg)  cos(cp)  +  [<I>o2/(47t2L)]  cp2  on  the  ground  band, 
and  kinetic  energy  q1 !2Cz-  Here,  <f>  is  the  magnetic  flux  in  the  loop.  Assuming  the  SSET 
stays  at  the  lowest  band,  where  I  ~  1 o  sin  (cp),  classical  dynamics  of  phase  in  the 
oscillator  is  thus  analogous  to  a  particle  moving  in  a  sinusoidally  modulated  parabolic 
potential.  The  dynamics  is  similar  to  that  of  the  rf-SQUID  [3],  i.e.,  a  single  junction 
shunted  with  a  loop  inductance,  except  that  since  our  loop  is  not  fully  superconducting, 
flux  quantization  or  flux  jumps  do  not  exist.  At  small  driving  amplitude  (linear  regime), 
the  phase  particle  experiences  harmonic  oscillations  around  cp  =  0  at  the  frequency /p  = 
1/(2ti)  (L  1 1  L]C)'112  which  is  controlled  by  gate-tuned  Lj.  This  mode  of  operation,  where 


the  L-SET  works  as  a  charge-to-firequency  transducer,  we  call  the  "plasma  oscillation 
mode". 


FIGURE  2.  Measured  frequency  response  of  the  reflection  coefficient  I  T  |  as  a  function  of  excitation 
power:  data  (left)  and  simulation  (right)  [4],  "Critical  power"  Pc  is  marked  by  the  arrow.  In  order  to 
guarantee  convergence  of  the  simulation,  the  calculations  were  performed  using  a  slightly  larger 
dissipation  than  present  in  the  experiment. 

At  higher  oscillation  amplitude,  in  the  “non-harmonic  regime”,  the  particle  sees  a 
different  local  curvature  of  the  potential  due  to  the  coscp  term,  thus  changing  the 
oscillation  frequency.  At  very  high  amplitude,  the  cosine  modulation  becomes 
effectively  averaged  out.  Thus,  when  increasing  excitation,  we  expect  a  change  of 
resonant  frequency  from  fvtofo-  1/(2ti)  (LC)'  .  This  change  of  resonance  frequency  at 
a  critical  excitation  power  Pc  is  seen  in  experimental  data  of  Fig.2.  The  dynamics  of  the 
L-SET,  however,  is  still  very  sensitive  to  the  %-induced  variations  of  the  cosine- 
potential  when  P  is  slightly  above  Pc- 

Charge  resolution  was  measured  using  both  amplitude  and  phase  read-out  without  any 
significant  difference  between  these  methods.  In  the  plasma  oscillation  mode  we  got  sq 
=  2.0-10'1 2 3e/v/Hz  at  the  maximum  power  Pc  — 116  dBm,  corresponding  to  20  fW 
dissipation  in  the  whole  resonator  circuit.  Note  that  due  to  Cooper  pair  tunneling,  the 
power  is  not  dissipated  in  the  SSET  island.  In  the  non-harmonic  mode,  significantly 
better  sensitivities  were  obtained:  sq  =  1.4-10"4e/VHz  was  measured  at  -100  dBm.  Due 
to  a  higher  impedance  mismatch,  power  dissipation  can  be  kept  below  20  fW  in  this 
case  as  well.  The  non-harmonic  mode  offers  good  possibilities  to  reach  uncoupled 
energy  sensitivity  of  e  ~  ft  using  an  optimized  sample  in  the  present  configuration  [4], 
Our  simulations  indicate  that  by  employing  a  sample  with  Rset  =  35  kO,  Ej=  0.45  K 
and  Lj  =  60  nH,  a  resolution  of  sq  ~  MO^e//  can  be  reached. 

Fruitful  discussions  with  T.  Heikkila,  G.  Johansson,  R.  Lindell,  H.  Seppa,  and  J.  Viljas 
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by  the  Large  Scale  Installation  Program  ULTI-3  of  the  European  Union. 
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We  present  an  experimental  characterization  in  the  quantum  regime  of  Nb/AlOx/Nb 
Josephson  devices.  Josephson  Junctions  and  folly  integrated  Josephson  devices  consisting  of  a  rf 
SQUID  coupled  to  a  readout  system  based  on  a  dc  SQUID  sensor  have  been  measured.  Data  on 
the  decay  rate  from  the  metastable  states  show  evidence  of  Macroscopic  Quantum  Tunneling1  with 
a  very  low  dissipation  level  and  a  high  insulation  of  the  probes  from  the  external  environment. 


Figure  1  Experimental  histograms 
of  the  switching  flux  distribution 
P((I>X)  at  different  temperatures.  The 
mean  value  of  the  flux  distributions 
is  reported  in  the  inset  as  a  function 
of  T,  together  with  the  theoretical 
prediction  in  thermal  (solid  line)  and 
quantum  limit  (dashed  line).  The  rf 
SQUID  parameters  are  L=80  pH, 
C=2.3  pF,  Ic=9  pA,  and  pL=2. 
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Figure  2.  Data  on  the  distribution  width 
a  as  function  of  T273  is  shown  for  the 
SQUID  system.  Experimental  data  in  the 
thermal  regime  are  compared  with  the 
theoretical  prediction  (dots  and  solid  line) 
obtained  for  R=1  MQ  and  C=2.3  pF.  We 
also  report  the  best  fit  line  of  the  data  in 
the  quantum  regime  (triangles  and  dashed 
line).  The  interception  of  this  line  with 
the  thermal  curve  gives  an  experimental 
extrapolation  of  the  crossover 
temperature:  Ts=250  mK. 


Furthermore  we  consider  the  theoretical  problem  of  resonant  phenomena  in  the  Macroscopic 
Quantum  Tunneling  for  an  rf-SQUID  in  presence  of  an  external  field1 2.  The  transition  probability 
between  energy  levels  is  studied  by  varying  the  parameters  of  the  potential  describing  the  system 
in  a  way  that  the  pumping  level  in  the  left  potential  well  is  close  to  some  level  in  the  right  potential 
well.  The  dependence  of  transition  probability  from  the  external  drive  of  the  system  shows  two 
resonance  peaks,  the  former  connected  with  resonance  tunneling  and  the  latter  with  resonance 
pumping.  Relative  position  of  peaks  depends  on  the  pumping  energy  and  on  the  potential 
parameters. 
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Figure  3.  Transition  probability  W  as  function  of  the  external  magnetic  flux  (px .  The  curve  is 
calculated  for  different  values  of  the  pumping  frequency.  The  two  peaks  shown  in  a)  and  b)  are 
connected  with  the  effects  of  resonance  tunneling  between  levels  in  different  wells  and  the 
resonance  pumping  between  levels  in  the  same  well.  As  predicted,  the  peaks  due  to  the  resonance 
pumping  moves  with  the  pumping  frequency  oo  until  it  superimposes  to  the  tunneling  peak,  as  it  is 
clear  in  Fig.  3b. 
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We  report  recent  experiments  that  demonstrate  quantum  tunneling  and  energy  level  quantization  of 
Josephson  vortices  at  temperatures  below  100  mK.  A  single  vortex  behaves  as  a  macroscopic  particle 
with  a  spatial  extent  of  several  micrometers  which  tunnels  through  a  potential  barrier  created  by  a 
magnetic  field  applied  to  the  annular  Josephson  junction  [1],  We  probe  the  quantum  properties  of  this 
particle  by  investigating  the  statistics  of  a  set  of  tunneling  events.  Measurements  of  vortex  escape  from 
a  potential  well  in  the  presence  of  microwave  radiation  clearly  indicate  discrete  energy  levels  of  the 
vortex.  Both  the  tunneling  rate  and  the  energy  level  separation  can  be  tuned  in  experiment  by  varying 
the  applied  magnetic  field. 

The  experiments  on  ultra-narrow  annular  long  junction  in  the  absence  of  initially  trapped  vortices  allow 
for  studying  the  properties  of  a  vortex  “molecule”,  consisting  of  a  bound  state  of  vortex  and  antivortex 
[2],  We  find  that  a  switching  from  the  superconducting  state  to  the  resistive  state  occurs  as  the 
dissociation  of  this  pinned  vortex-antivortex  pair.  The  interaction  potential  of  the  bound  state  is 
controlled  by  external  magnetic  field  and  dc  bias  current.  The  figure  below  shows  our  experimental 

data  (dots)  for  the  lowest  energy  level 
separation  of  vortex-antivortex  bound 
state  measured  by  microwave 
spectroscopy  within  10-20  GHz 
frequency  range.  Solid  line  is  the 
calculated  vortex  “molecule”  level 
separation  [2];  dashed  line  is  the 
standard  prediction  of  the  lumped 
junction  model.  The  inset  schematically 
shows  the  sample  geometry. 

Engineering  an  energy  profile  for  a 
vortex  in  a  long  Josephson  junction 
opens  an  opportunity  of  designing 
vortex  qubits.  The  energy  profile  for  the 
vortex  can  be  created,  e.g.,  by  making 
the  junction  of  a  particular  geometrical 
shape  and  applying  an  external  magnetic 
field  [3].  Possible  coupling  schemes  of 
vortex  qubits  [4]  are  feasible  for  creating  and  monitoring  entanglement  in  multi-qubit  circuits.  Vortex 
qubits  allow  for  a  convenient  state  preparation  and  readout  using  superconducting  Rapid  Single  Flux 
Quantum  (RSFQ)  circuitry  [5], 
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We  have  investigated  the  statistical  behavior  of  Josephson  junctions  with  respect 
to  their  switching  from  the  zero-voltage  state  when  ac-perturbations  are  applied.  The 
investigations  are  inspired  by  the  experimental  observations  of  anomalous  switching 
distributions  reported  in  for  temperatures  below  the  quantum  transition  temperature 
[1,  2].  These  measurements  have  produced  multi-peaked  switching  distributions  when 
microwaves  are  applied  to  the  system,  and  the  results  have  been  interpreted  as  a 
signature  of  intrinsic  quantum  energy  levels  excited  by  the  external  microwaves.  Our 
work  has  been  focused  on  the  classical  limit  of  the  system,  where  thermal  fluctuations 
dominate  quantum  uncertainties,  and  where  the  governing  model  equations  of  motion 
can  be  written  as  stochastic  differential  equations. 

Extensive  simulation  results  reveal  [4,  5]  that  the  entirely  classical  models  can 
reproduce  the  multi-peaked  switching  distribution  signatures  observed  experimentally 
in  the  low  temperature  limit.  We  have  investigated  both  the  small  area  Josephson 
junction  and  the  magnetic  field  embedded  long  annular  Josephson  junction  with  a 
single  trapped  magnetic  fluxon.  The  simulated  current- voltage  characteristics  of  the 
latter  system  can  with  good  approximation  be  obtained  from  a  simplified  effective 
pendulum  equation,  and  the  characteristics  of  the  small  and  long  annular  junctions 
are  therefore  very  similar  [3]. 

The  essence  of  the  system  is  summarized  in  Figure  1  (left),  where  <f>  is  the  quantum 
mechanical  phase-difference  between  the  superconductors  of  the  junction.  The  escape 
is  assisted  by  two  components:  1)  the  thermal  bath,  and  2)  forced  oscillations  from 
the  ac-perturbation.  We  exemplify  our  investigations  of  the  small  area  Josephson 
junction,  modeled  by  a  pendulum  equation,  through  a  set  of  switching  distributions 
shown  in  Fig.  1.  The  middle  plot  shows  the  switching  statistics  for  several  different 
applied  frequencies.  The  locations  of  the  secondary  (left)  probability  peak  as  a  func¬ 
tion  of  driving  frequency,  fl,  is  shown  on  top  along  with  the  results  of  harmonic  and 
anharmonic  analysis  of  the  intrinsic  resonance  as  a  function  of  bias  current,  rj.  As  is 
well  known,  the  linear  resonance  frequency  of  the  system  is  a  function  of  the  bias, 
rj.  Thus,  as  the  bias  is  swept  (slowly)  from  low  values  toward  the  critical  value,  the 
resonance  frequency  of  the  system  is  tuned  from  the  plasma  frequency  to  zero.  The 
classical  interpretation  of  the  anomalous  switching  statistics  is  that  if  the  driving  fre¬ 
quency  (or  a  significant  sub  or  super  harmonic)  is  less  than  the  plasma  frequency,  then 
multi-peaked  switching  distributions  may  arise  by  simple  resonant  coupling  between 
the  inherent  system  resonance  and  the  external  frequency,  provided  the  ac-amplitude 
is  tuned  appropriately. 

Similarly  to  the  small  area  junction,  we  can  demonstrate  the  same  behavior  for 
long  annular  junctions  in  a  magnetic  field  and  a  single  trapped  fluxon.  The  system  is 
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here  modeled  by  the  sine-Gordon  equation.  Figure  1  (right)  shows  anomalous  switch¬ 
ing  distributions  for  the  escape  of  the  trapped  fluxon  in  the  magnetic  field  induced 
potential  well.  Again  we  observe  anomalous  multi-peaked  switching  distributions, 
generated  by  the  ac-perturbation,  and  we  verify  the  interpretation  of  the  anomalous 
distributions  as  a  direct  result  of  a  simple  forced  oscillation,  which  resonates  near 
specific  values  of  the  bias  current. 


Figure  1:  (left)  Sketch  of  the  energetics  of  the  system  under  consideration:  an  object 
is  escaping  a  potential  well  due  to  a  combination  of  forced  oscillations  and  thermal 
fluctuations,  (middle  -  see  [j]  for  details)  Escape  distributions  for  different  frequen¬ 
cies  u,  applied  to  a  pendulum  model  of  a  small  area  Josephson  junction,  as  a  function 
of  bias  current,  rj.  Upper  part  shows  the  agreement  with  harmonic  (solid)  and  anhar- 
monic  (dashed)  analyses,  (right  -  see  [5]  for  details)  Escape  distributions  for  different 
frequencies  u>,  applied  to  a  sine-Gordon  model  of  a  long  annular  Josephson  junction 
with  a  single  trapped  fluxon,  as  a  function  of  bias  current,  rj.  Upper  part  shows  the 
agreement  with  harmonic  (solid)  and  anharmonic  (dashed)  analyses. 

We  will  present  comprehensive  numerical  simulations  of  the  classical  system,  an¬ 
alytical  perturbation  results  of  the  nonlinear  equations  to  substantiate  the  data,  and 
some  comparisons  to  experimental  measurements  conducted  at  temperatures  well 
above  the  quantum  transition. 
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Recently,  different  types  of  qubits,  all  based  on  Josephson  junctions,  have  been 
experimentally  demonstrated.  Flux  [1,2],  phase  [3,4,5]  and  phase-charge  [6]  qubits  have 
been  operated  as  single  devices,  while  charge  qubits  have  also  been  used  in  an  entangled 
couple,  showing  quantum-coherent  behaviour  [7]  and  operation  as  conditional  gate  [8]. 
Which  one  of  the  proposed  qubits  will  be  the  most  promising  in  view  of  the  implementation 
in  a  practical  quantum  computer  is  still  an  open  issue.  The  main  points  that  concerns  a 
feasible  qubit  system  relay  on  the  following  features:  the  single  qubit  should  have  a  long 
decoherence  time,  its  writing  and  reading  system  must  be  simple,  fast  and  should  allow  a 
one  shot  detectability  with  high  efficiency.  The  single  qubit,  moreover,  should  allow  a 
simple  and  switchable  connection  with  one  or  more  qubits  to  realize  the  core  of  a  C-NOT 
gate. 

In  our  talk  we  will  present  the  system  we  have  developed  to  realize  single  and  coupled  set 
of  qubits.  At  present  we  have  two  systems  that  may  operate  as  a  phase  and  flux  qubits.  The 
flux  qubit  is  an  rf  SQUID  realized  with  standard  Niobium  trilayer  technology,  read  by  a 
hysteretic  dc  SQUID  working  as  a  one  shot  detector  of  the  rf  SQUID  flux  states.  The 
hysteretic  dc  SQUID  may  operate  also  as  a  phase  qubits,  driven  by  short  pulses  of 
microwaves.  On  this  system  we  made  a  set  of  measurements  to  show  induced  oscillations 
between  the  first  two  energy  levels  of  the  dc  SQUID. 

A  third  scheme  to  implement  a  flux  qubit,  where  the  reading  system  is  integrated  with  the 
qubit  has  been  realized  and  tested.  We  will  present  the  main  idea  underlying  the  qubit 
design  and  some  preliminary  test  on  the  system,  showing  the  correct  behaviour  of  the 
proposed  device. 

To  realize  an  efficient  coupling  between  flux  qubits  we  have  proposed  and  realized  a 
switchable  coupling  making  use  of  the  variable  inductance  of  a  hysteretic  dc  SQUID  driven 
by  short  pulses  of  current.  This  coupling  transformer  looks,  at  the  moment,  one  of  the  most 
efficient  schemes  to  implement  a  flux  based  c-not  gate. 
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Superconducting  devices  based  on  the  Josephson  effect  have  proved  to  be  promising 
candidates  for  the  implementation  of  quantum  computing.  After  a  demonstrative  phase, 
concerning  the  capability  to  implement  and  couple  Josephson  qubits  [1-8],  the  actual 
efforts  concern  the  design  and  realization  of  new  optimised  elements  with  small 
decoherence,  efficient  readout,  simple  and  complete  coherent  manipulation,  integration 
and  coupling  of  more  qubits.  In  this  direction,  we  introduce  a  flux  qubit  based  on  a  fully 
gradiometric  double  rf  SQUID,  with  a  readout  scheme  based  on  the  direct  injection  of 
current  in  a  large  junction  in  the  SQUID  loop.  The  device  has  been  realized,  and  we 
present  a  preliminary  experimental  characterization  at  4.2  K  showing  the  capability  to 
perform  an  efficient  flux  state  readout. 

The  qubit  is  essentially  an  rf  SQUID  with  the  superconducting  loop  interrupted  by  a 
very  small  inductance  hysteretic  dc  SQUID  [9],  The  system  is  described  by  the 
magnetic  flux  threading  the  large  loop,  with  a  potential  shaped  as  a  double  well  (if  the 
externally  applied  flux  is  near  a  half  flux  quantum,  and  for  appropriate  constructive 
parameters).  The  distinct  flux  states  in  the  two  minima  can  be  used  as  computational 
states  “0”,  “1”.  The  flux  applied  to  the  large  loop  controls  the  potential  symmetry,  and 
can  be  used  to  prepare  the  system  in  one  of  the  two  flux  states,  and  to  manipulate  the 
relative  phase  between  them  during  the  computation.  The  flux  applied  to  the  small  loop 
controls  the  potential  barrier  height,  allowing  the  freezing  or  the  “rotation”  between  the 
flux  states  [10],  The  manipulation  can  be  performed  by  means  of  non-adiabatic 
changes,  realized  for  example  by  using  Rapid  Single  Flux  Quantum  Logic  elements  [11, 
12].  Our  device  is  designed  in  order  to  have  a  gradiometric  configuration  for  both  the 
small  and  the  large  loops.  This  has  the  double  advantage  to  protect  the  qubit  against 
magnetic  noise,  and  to  allow  the  device  self-bias:  if  a  single  flux  quantum  (or  an  odd 
number  of  flux  quanta)  is  trapped  in  the  gradiometric  large  loop,  the  qubit  is 
automatically  biased  near  the  symmetric  condition,  required  for  the  computation.  The 
readout  scheme  is  similar  to  that  used  for  the  Quantronium  charge-phase  qubit  [4],  and 
consists  in  a  large  junction  inserted  in  the  qubit  loop.  A  current  pulse  directly  applied  to 
this  junction  can  induce  a  transition  to  the  voltage  state  for  a  critical  value  that  strongly 
depends  on  the  qubit  flux  state.  This  allows  an  effective  one-shot  qubit  state  readout. 


In  the  above  figure  are  shown:  (a)  the  simplified  scheme  of  the  considered  flux  qubit 
(for  simplicity  it  is  not  considered  the  gradiometric  configuration);  (b)  the  plot  of  the 
readout  critical  current  Ic  for  different  fluxes  &x  applied  to  the  large  loop  (horizontal 
axis),  measured  at  4.2  K.  The  different  curves  correspond  to  different  fluxes  applied  to 
the  small  control  loop  (0c),  related  to  different  barrier  height  and  distances  between 
minima.  Since  the  measured  spread  of  the  critical  current  due  to  thermal  noise  is  below 
1  pA,  the  one-shot  readout  can  be  performed  also  at  4.2  K  with  a  very  high  efficiency. 
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We  show  that  in  general  the  transition  from  the  classical  to  quantum  behavior  depends 
on  the  probing  length  scale,and  occurs  for  microscopic  length  scales,  except  when  the 
interactions  between  the  particles  are  very  weak.  This  transition  explains  why,  on 
macroscopic  length  scales,  physics  is  described  by  classical  equations. 

In  particular,  the  smooth  transition  from  localized  quantum  mechanical  information  to 
classical  waves  may  allow  one  to  holographically  reconstruct  entangled  quantum 
states. 
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We  examine  the  effect  of  multilevels  on  decoherence  and  dephasing  properties  of  a 
quantum  system  consisting  of  a  non-ideal  two  level  subspace,  identified  as  the  qubit  and  a 
finite  set  of  higher  energy  levels  above  this  qubit  subspace.  The  whole  system  is  under 
interaction  with  an  environmental  bath  through  a  Caldeira-Leggett  type  coupling.  The 
model  that  we  use  is  an  rf-SQUID  in  the  macroscopic  quantum  coherence  regime  and 
coupled  inductively  to  a  flux  noise  characterized  by  an  environmental  spectrum.  The  model 
interaction  can  generate  dipole  couplings  between  the  qubit  subspace  and  a  high  number  of 
excitation  levels.  The  decoherence  properties  of  the  qubit  subspace  is  examined 
numerically  using  the  master  equation  formalism  of  the  system's  reduced  density  matrix. 
We  numerically  examine  the  relaxation  and  dephasing  times  as  the  environmental 
frequency  spectrum,  the  environmental  temperature,  and  the  multilevel  system  parameters 
are  varied.  We  observe  that,  these  time  scales  receive  contribution  from  all  available 
energies  in  the  noise  spectrum  (even  well  above  the  system's  energy  scales)  stressing  the 
dominant  role  played  by  the  non-resonant  (virtual)  transitions.  The  relaxation  and 
dephasing  times  calculated,  strongly  depend  on  the  number  of  active  levels.  Under  the 
influence  of  these  effects,  the  validity  of  the  two  level  approximation  is  restricted  not  by  the 
temperature  but  by  the  dipole  couplings  as  well  as  the  availability  of  the  environmental 
modes  at  low  temperatures. 

To  stress  the  importance  of  the  non-resonant  transitions,  we  use  other  model 
spectral  profiles  and  study  the  effect  of  the  central  frequency,  the  width  and  the  total 
spectral  area  on  the  decoherence  times.  We  find  that,  due  to  the  strong  influence  of  the  non¬ 
resonant  transitions,  the  spectral  parameters  affect  the  decoherence  times  only  through  the 
dependence  on  the  spectral  area.  In  this  regard  the  weight  of  the  resonant  contributions  is 
negligible  as  compared  to  that  of  the  entire  noise  spectrum. 

The  credibility  of  our  results  depends  on  the  accuracy  of  the  Bom  approximation  in 
the  short  time  evolution.  In  this  regard,  it  is  crucial  that  the  approximation  holds  in  the 
range  of  the  numerical  data.  This  has  been  also  confirmed  by  checking  the  renormalization 
of  the  noise  spectrum  at  the  RPA  level. 


Figure  1.  Relaxation  and  dephasing  rates 
against  the  number  of  levels  for  different 
spectral  widths  at  T=0.  The  open  and  solid 
symbols  refer  to  dephasing  and  relaxation 
times  respectively. 


Figure  2.  Decoherence  times  for  rectangular 
and  Lorentzian  spectra 


Figure  3.  Renormalization  of  the  noise 
spectrum  at  the  RPA  level 
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The  superconductor-insulator  (SI)  transition  in  a  single  Josephson  junction  (JJ)  and  JJ 
arrays  in  one  and  two  dimensional  geometry  have  been  extensively  studied.  Experimentally 
observed  phase  diagrams  in  the  E/Ec  -  g  plane  show  two  features  [1-5],  Here  Ej,  Ec  and 
g=Ro/Rs  are,  respectively,  the  Josephson  energy,  capacitive  Coulomb  energy  of  the 
superconducting  islands,  and  the  Ohmic  shunt  resistance  in  units  of  the  quantum  resistance, 
RQ=h/4e2.  One  is  dissipative  transition  controlled  by  g  and  the  other  is  quantum  phase 
transition  controlled  by  Ej/Ec.  Here  our  focus  is  on  the  Ej/Ec  axis  at  g=0  where  the  critical 
Ej/Ec~3  for  a  single  JJ  and  ~0.3  for  the  48x40  JJ  arrays.  This  order  of  magnitude  difference 
was  argued  to  be  an  evidence  of  cooperative  phenomena  and  the  phase  boundary  in  the 
48x40  JJ  arrays  a  true  phase  boundary. 

The  theoretically  irrefutable  treatment  of  the  fully  2D  quantum  mechanical  and  strongly 
correlated  systems  near  the  phase  boundary  is  an  ever-lasing  challenge.  One  simplification 
is  to  neglect  the  quantum  fluctuation  due  to  the  capacitive  Coulomb  energy  and  assume  that 
the  single  parameter  in  the  resulting  2D  XY  model  Ej/T  where  T  is  temperature  may  play 
the  role  of  E/Ec  at  T=0.  The  difference  between  the  fully  quantum  mechanical  JJ  and  the 
classical  2D  XY  model  is  that  the  former  undergoes  a  quantum  phase  transition  (QPT) 
whereas  the  latter  a  classical  phase  transition  (CPT).  The  qualitative  resemblance 
nevertheless  lies  in  the  fact  that  both  undergo  the  Berezinskii-Kosterlitz-Thouless  (BKT) 
transition.  More  importantly,  the  strong  size  dependence  seems  to  be  rather  universal  in  two 
dimension.  In  fact,  while  the  2D  XY  model  has  been  studied  extensively  in  the  past,  there 
are  two  noteworthy  development  recently  concerning  the  strong  size  dependence.  One  is  in 
the  context  of  the  planar  magnets.  It  is  well-known  that  in  spite  of  the  fact  that  the  BKT 
transition  describing  the  SI  transition  is  infinite  order,  in  a  dozen  of  planar  magnets,  the 
magnetization  appears  like  in  the  2-nd  order  transition  with  a  universal  magnetization 
exponent  -0.23.  Bramwell  and  Holdsworth  [6]  used  renormalization  group,  Monte  Carlo 
and  finite-size  scaling  to  show  that  the  transition  looks  like  2-nd  order  with  Tc~1.08  for  103 
spins  and  -1.02  for  104  spins,  which  are  fairly  modified  from  Tc=0.892  believed  to  be  exact 
for  the  infinite  system.  The  size  dependence  is  significant,  reflecting  the  peculiar 
singularity  of  the  correlation  length  in  the  BKT  transition.  The  other  development  is  the 
transfer-matrix  density-matrix-renormalization-group  (TM-DMRG)  method  by  the  present 
author  [7],  Using  the  Roomany-Wyld  finite-size  scaling,  without  any  unknown  parameters, 
nor  relying  on  the  BKT-RG,  we  found  that  Tc=l  .07  for  2000-3000  spins  in  agreement  with 
[6] 

To  account  for  the  strong  size  dependence  in  the  critical  Ej/Ec~Ej/T,  we  need  to  go  one 
step  further  and  look  for  a  signature  of  the  SI  transition  as  a  function  of  the  system  size. 
Note  that  the  2D  XY  model  is,  under  Villain  approximation,  equivalent  to  the  ID  quantum 
sine-Gordon  (SG)  model,  with  the  correspondence  Ej/Ec  -  1/p2  [8], 


1 


"  =  £{• 

t=1 


The  ground  state  of  the  SG  Hamiltonian  can  be  analyzed  by  DMRG.  We  have  found  that 
there  is  a  spontaneous  symmetry  breaking  (SSB)  of  the  ground  state  as  a  function  of  the 
system  size  [9].  The  identification  of  the  BKT  transition  as  a  SSB  is  then  demonstrated  by 
combining  DMRG  with  the  Roomany-Wyld  finite-size  scaling.  Fig.l  is  the  ground  state 
phase  diagram  in  the  J32-L  plane.  With  increase  of  the  size,  Ej/Ec  «  1/p2  decreases  like 
p2(L=5)/p2(L=50)  -0.1  which  is  consistent  with  the  experimental  findings. 

After  this  is  done,  it  is  now  clear  that  we  can  accurately  analyze  the  size  dependences  of 
(1)  the  QPT  in  the  JJ  chain  with  on-site  Coulomb  fluctuation  and  (2)  the  CPT  in  the  2D  XY 
model  for  the  JJ  arrays  with  relatively  large  superconducting  islands.  We  will  report  on  the 
size  dependence  of  the  superfluid  densities  of  these  systems. 


FIG.l  The  ground-state  phase  diagram  in  the  p2-L  plane 
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Isolated  large  area  Josephson  junctions  can  be  used  as  qubits  for  quantum  computation 
[1-3].  Readout  of  the  junction  can  be  performed  by  allowing  the  system  to  tunnel  to  the 
voltage  state.  Single  qubit  gate  operations  are  performed  by  capacitively  coupling  short 
microwave  pulses  into  the  junction  bias  current.  These  pulses  must  be  carefully 
designed  so  as  to  achieve  purely  ox  or  oy  operations.  The  large  capacitance  of  the 
Josephson  junctions,  typically  a  few  pF,  allows  for  the  use  of  capacitive  inter-qubit 
coupling  [4,5].  While  fixed  physical  capacitors  can  be  used  to  couple  junctions,  a 
tunable  capacitance  coupling  scheme  not  only  simplifies  gates,  but  also  reduces 
perturbation  of  one  qubit  by  the  tunneling  of  the  other  qubit  to  the  voltage  state.  In  this 
talk,  I  will  discuss  challenges  in  implementing  single  and  two  qubit  gates  and  a  tunable 
capacitor  coupling  scheme  based  on  a  Cooper  pair  box. 

The  microwave  pulses  for  single  qubit  gate  operations  can  be  generated  by  room 
temperature  electronics.  Ideal  single  qubit  gate  pulses  have  a  roughly  Gaussian 
envelope  and  fixed  phase  and  have  durations  of  several  nanoseconds  [6].  To  generate 
these  pulses  one  needs  full  vector  control  of  the  microwave  source,  that  is  control  in 
both  amplitude  and  phase.  Errors  introduced  in  the  pulse  phase  by  the  amplitude 
modulation,  as  typically  occurs  in  most  pulse  modulators,  results  in  single  gate 
operations  that  can  be  described  as  sequenced  application  of  ox  then  cry  then  crx  on  the 
qubit.  Without  proper  control,  this  results  in  effectively  a  reduced  modulation 
amplitude  during  Rabi  oscillations.  Also,  since  the  microwave  drive  pulses  are  in  the 
strong  drive  regime,  the  rotating  wave  approximation  no  longer  applies  [7]  and  the 
envelope  of  the  microwave  pulse  must  be  synchronized  to  the  microwave  phase  with 
jitter  of  less  than  10  ps  for  typical  qubit  parameters  and  desired  gate  fidelities. 
Generating  these  desired  microwave  pulses  is  possible  using  state  of  the  art  electronics. 

In  the  case  of  fixed  capacitive  coupling,  two  qubit  operations  can  be  performed  by  non- 
adiabatically  biasing  the  qubits  to  tune  the  energy  level  spacing  of  one  qubit  into 
resonance  with  a  second  qubit.  These  bias  manipulation  pulses  must  maintain  the  high 
signal  to  noise  ratio  of  the  bias  current  while  simultaneously  having  a  high  bandwidth. 
After  manipulation,  readout  of  one  of  the  qubits  by  allowing  it  to  tunnel  to  the  voltage 
state  results  in  a  burst  of  current  to  any  qubit  coupled  to  it.  I  will  briefly  describe  an 
experiment  that  measures  the  resulting  correlated  switching  between  the  two  qubits  for 
the  case  of  fixed  capacitive  coupling  This  effect  can  be  reduced  by  using  a  tunable 
capacitor  between  the  qubits  that  can  be  turned  off  during  measurement. 

A  tunable  capacitor  can  be  implemented  by  using  a  Cooper  pair  box  (CPB).  The 
capacitance  of  a  typical  CPB  for  this  application  can  be  tuned  from  0  fF  to  roughly  50 
fF  through  gate  voltage  control.  For  the  simplest  version  of  the  tunable  capacitor  to 
work,  one  must  avoid  resonant  interactions  between  the  qubit  and  the  CPB  so  as  to  keep 


the  CPB  in  the  ground  state.  This  can  be  achieved  by  designing  the  Josephson  coupling 
energy  Ej,  of  the  CPB  to  be  higher  than  the  energy  spacing  of  the  qubit,  or  by  only 
limiting  die  gate  voltage  to  the  regime  where  the  CPB  energy  spacing  is  more  than  the 
qubit  energy  spacing.  One  must  also  keep  the  charging  energy  Ec,  larger  than  Ej  to 
maintain  a  large  on/off  ratio  of  capacitance.  This  tunable  capacitor  scheme  can 
potentially  be  combined  with  a  bus  bar  to  allow  coupling  between  non-nearest 
neighbors. 
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Experimental  Studies  of  a  Josephson  Charge  Qubit  with  RF-SET  read-out. 
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41296  Gdteborg,  Sweden 

In  order  to  experimentally  realize  a  superconducting  charge  qubit,  one  must  overcome  several 
challenging  problems.  The  first,  which  is  surprisingly  subtle  and  difficult,  involves  producing 
Cooper-pair-boxes  with  robust  even-parity  states.  One  must  then  optimize  the  readout  while 
simultaneously  minimizing  its  back-action  on  the  qubit.  Finally,  the  mechanisms  of  decoherence 
in  Josephson  junction  qubits  must  be  understood  if  we  are  to  scale  up  to  many  qubits.  We  must 
identify  both  the  sources  of  dephasing  and  relaxation,  and  how  they  couple  to  the  qubit.  Although 
there  is  much  to  learn  from  the  mature  field  of  spin  physics,  superconducting  qubits  are 
fundamentally  different  in  that  they  are  mesoscopic  objects,  the  noise  is  very  asymmetric,  and 
measurements  are  not  made  on  ensembles. 

We  present  experimental  studies  of  a  Josephson  charge  qubit  (Cooper-pair-box)  that  uses  radio- 
frequency,  single-electron-transistor  (RF-SET)  readout.  We  describe  the  evolution  of  our  sample 
fabrication  process  that  led  to  robust  even-parity  boxes  with  SET's  that  can  be  operated  in  the 
sub-gap  regime  in  a  manner  that  minimizes  the  backaction.  Fast  DC  pulses  are  used  to 
manipulate  the  qubit  and  coherent  oscillations  with  a  high  fidelity  are  consistently  observed.  The 
reduction  from  100%  is  understood  as  an  adiabatic  effect  of  the  finite  pulse  rise-time. 
Measurements  of  dephasing  for  the  full  range  of  gate  charge  and  on  time  scales  extending  to  the 
sub-nanosecond  regime  have  been  made.  These  allow  us  to  attribute  the  main  source  of 
dephasing  to  fluctuations  that  couple  predominantly  to  charge. 

Measured  relaxation  times  vary  considerably  and  can  be  smaller  than  what  is  expected  from 
either  the  Ohmic  control  circuitry  or  extrapolation  of  the  known  1/f  spectrum  to  GHz  frequencies. 
Investigations  of  various  relaxation  processes  as  well  as  backaction  due  to  the  SET  are  also 
presented. 
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Josephson  junction  qubit  network  with  current-controlled  interaction 
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After  several  years  of  great  experimental  achievements  on  Josephson  junction  (JJ) 
qubits  [2-4],  coupling  of  qubits  is  now  becoming  feasible.  A  few  challenging 
experiments  with  coupled  JJ-qubits  have  been  reported  [5].  However,  so  far 
experiments  on  coupled  JJ  qubits  have  been  performed  without  direct  physical  control 
of  the  qubit-qubit  coupling. 

We  present  a  new  solution  for  controllable  physical  qubit-qubit  coupling  [1].  The 
network  has  the  following  properties:  (a)  nearest-neighbour  qubit-qubit  coupling 
controlled  by  external  bias  current,  (b)  qubits  parked  at  the  degeneracy  points,  also 
during  qubit-qubit  interaction,  (c)  separate  knobs  for  controlling  individual  qubits  and 
qubit-qubit  coupling,  (d)  scalability.  An  important  feature  is  that  the  network  is  easily 
fabricated,  and  is  in  line  with  current  mainstream  experiments. 


Network  of  loop-shaped  SCPT  charge  qubits,  coupled  by  large  Josephson  junctions. 
The  interaction  of  the  qubits  (i)  and  (i+1)  is  controlled  by  the  current  bias  hi  or  by 
simultaneous  current  biasing  of  readout  junctions.  Individual  qubits  are  controlled  by 
voltage  gates,  ngi=CgiVgi.  The  Josephson  energies  of  the  coupling  JJs,  E] ,  and  the 
readout  JJs,  Ej ,  are  much  larger  than  the  charging  energies,  Ej(m)  »  E^m) . 

The  network  under  consideration  consists  of  a  chain  of  charge  qubits  -  Single  Cooper 
Pair  Transistors  (SCPT)  -  with  loop-shaped  electrodes  coupled  together  by  current 
biased  coupling  JJs  at  the  loop  intersections.  The  loop  design  creates  an  (inductive) 
interface  to  the  qubit  by  means  of  circulating  currents  [4].  We  employ  these  current 
states  in  the  qubit  loops  to  create  controllable  coupling  of  neighbouring  qubits.  Left 
without  any  external  current  biasing  of  the  coupling  and  readout  JJs,  the  network  acts 
as  a  quantum  memory  of  independent  qubits  (neglecting  a  weak  residual  interaction 
[1]).  When  a  bias  current  is  sent  through  the  coupling  JJ,  the  current-current 
interaction  between  the  neighbouring  qubits  is  switched  on  and  increases  with 
increasing  bias  current.  Moreover,  if  both  of  the  readout  JJs  of  the  same  qubits  are 
biased  well  below  threshold,  again  there  is  nearest-neighbour  coupling  via  the 
circulating  currents.  However,  if  one  of  the  readout  JJs  is  current  biased,  this  only 
affects  that  particular  qubit  and  allows  readout  of  individual  qubits.  Thus  the  bias 
currents  serve  as  the  interaction  control  knobs.  The  loop  inductances  are  assumed  to  be 


sufficiently  small  to  neglect  qubit-qubit  coupling  via  induced  magnetic  flux,  as  well  as 
undesirable  qubit  coupling  to  the  magnetic  environment.  In  addition,  we  assume 
negligible  capacitive  coupling  between  the  islands,  which  are  well  shielded  by  the 
injection  leads. 

The  qubit  network  dynamics  in  the  charge  qubit  regime  Ec»Ej  is  described  by  the 
Hamiltonian 

x  £  _ 

H  =  Xhf-O  - ng> K;  - EJi  cos^o-,,  +  rii(Txi<7x{M)\  +  £*><7*0^, 

*=i  z  i*} 


which  consists  of  the  individual  qubits,  where  Oj  is  the  mean  phase  difference  induced 

by  current  bias  in  the  coupling  and  readout  JJs,  the  controllable  interaction,  and  finally 
a  weak  residual  coupling.  The  residual  coupling  effectively  connects  all  the  qubits,  but 
is  a  factor  hcop  IE)  « 1 ,  where  &p  is  the  coupling  JJ  plasma  frequency,  smaller  than 
the  maximum  controllable  coupling. 

In  order  to  couple  exclusively  the  qubits  (i)  and  (i+1)  one  should  apply  a  current  bias 
/*  ,  while  l)±t  -  0 .  In  this  case  the  qubit  coupling  amplitude  is  given  by  the  equation 

ft 


7, 


EjjEj(M) 
4  E)  cos^/' 


sin 


The  coupling  is  quadratic, 


(/*  / 1)  J ,  for  small  currents  and  diverges  when  the 


current  approaches  the  critical  current  Ibc .  An  alternative  way  to  switch  on  the  qubit 

interaction  is  to  apply  current  bias  (below  the  critical  current)  simultaneously  to  two 
neighbouring  readout  JJs,  which  will  give  a  similar  but  weaker  interaction. 

By  applying  current  bias  pulses  with  controlled  amplitude  and  duration  it  is  possible  to 
perform  quantum  gates  on  the  qubit  network,  eg.  control-NOT  and  cNOT-SWAP 
gates.  Readout  can  be  performed  either  using  the  readout  JJ:s  [6]  or  by  measuring  the 
charge  on  the  qubit  islands  [7]. 

Finally  we  emphasize  that  the  present  analysis  and  design  is  equally  relevant  in  the 
region  of  Ec~Ej,  characterizing  the  ’’Quantronium”  [4]. 

In  conclusion,  the  present  scheme  provides  a  realistic  solution  for  easy  local  control  of 
the  physical  coupling  of  charge  qubits  via  current  biasing  of  coupling  Josephson 
junctions  or,  alternatively,  pairs  of  readout  junctions.  The  design  is  in  line  with 
experimental  mainstream  development  of  charge  qubit  circuits  and  can  easily  be 
fabricated  and  tested  experimentally.  The  tunable  coupling  of  the  qubit  chain  allows 
easy  implementation  of  CNOT  and  CNOT-SWAP  operations. 
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Continuous  Monitoring  of  Rabi  Oscillations  in  a  Josephson  Flux  Qubit 
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Under  resonant  irradiation,  a  quantum  system  can  undergo  Rabi  oscillations  in  time,  a 
signature  of  quantum  coherence.  We  have  seen  evidence  for  such  oscillations  [1]  in  a 
continuously  observed  three-Josephson-junction  flux  qubit  [2] — a  macroscopic, 
SQUED-type  device  consisting  of  a  small  superconducting  ring  including  the  three 
junctions,  in  which  the  states  correspond  to  the  two  possible  directions  of  circulating 
supercurrent. 

In  any  discussion  of  quantum  phenomena,  prominent  attention  must  be  given  to 
the  readout  method.  Our  qubit  is  coupled  to  a  high-quality  (0t  ~  2000)  superconducting 
LC  oscillator  (i.e.,  a  tank  circuit)  tuned  to  the  Rabi  frequency,  acting  as  a  weakly 
coupled  but  sensitive  antenna  for  the  oscillations  in  magnetic  flux;  see  Fig.  1  for  the 
circuit  diagram,  and  Fig.  2  for  a  micrograph.  Since  the  Rabi  frequency  is  at  least  two 
orders  of  magnitude  smaller  than  the  one  of  the  external  field  (in  resonance  with  the 
qubit  level  splitting),  the  latter  is  effectively  decoupled  from  the  readout  circuit.  In 
addition  to  simplicity,  this  method  of  Rabi  spectroscopy  enables  a  long  coherence  time 
of  about  2.5  ps  (corresponding  to  an  effective  qubit  quality  factor  QRahl  ~  7000),  a 
record  for  this  class  of  devices.  In  conclusion,  our  recent  experiments  on  two  coupled 
flux  qubits  [3]  will  be  briefly  discussed. 
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Fig.  1:  Measurement  setup.  The  flux  qubit  is  inductively  coupled  to  a  tank  circuit.  The 
DC  source  applies  a  constant  flux  <I>e  «  <f>0/2.  The  HF  generator  drives  the  qubit  through 
a  separate  coil  at  a  frequency  close  to  the  level  separation  A/h  =  868  MHz.  The  output 
voltage  at  the  resonant  frequency  of  the  tank  is  measured  as  a  function  of  HF  power. 


Fig.  2:  The  aluminum  qubit  inside  the  niobium  pancake  coil. 
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We  investigate  the  dynamics  of  a  scalable  quantum  computer  employing  Josephson 
junction  (JJ)  prism  qubits  coupled  to  resonant  LC  buses.  The  basic  qubit  states  are  the 
symmetric  ground  |0)  and  antisymmetric  excited  |1)  states  of  a  two  triangular  cell  JJ 
prism  qubit.  The  physical  qubits  are  paired  into  logical  qubits  that  operate  in  a 
decoherence  free  subspace.  The  circulating  currents,  corresponding  to  the  |0)  and  |1) 
physical  qubit  states,  follow  the  perimeter  of  the  qubit,  and  are  symmetric  with  respect 
to  the  flux  threading  the  two  cells  of  the  qubit.  A  coupled  SQUID  gradiometer  that  is 
antisymmetric  with  respect  to  the  flux  threading  the  two  cells  will  not  register  any 
coupled  flux  until  a  Hadamard  gate  is  carried  out  on  the  physical  qubit,  immediately 
prior  to  a  measurement.  The  two  basic  states  are  then  transformed  into  ±  combination 
states  that  have  antisymmetric  ‘figure  8’  circulating  currents  whose  flux  is  then 
measurable  by  the  SQUID.  The  coupling  of  each  qubit  to  an  LC  resonant  bus  is  made 
by  means  of  a  loop  that  encircles  both  cells  of  the  qubit  symmetrically.  This  prevents 
the  direct  mutual  inductance  coupling  of  the  bus  to  both  antisymmetric  microwave  gate 
pulses  and  the  measuring  SQUID  gradiometer. 

One  and  two  logical-qubit  gates  and  bus  transfers  are  carried  out  using  bichromatic 
Molmer-Sorensen[l]  gates.  Since  a  resonant  LC  bus  can  only  accommodate  on  the 
order  of  50-100  logical  qubits,  the  system  can  be  enlarged  by  linking  many  buses 
together  in  an  open  network.  Two  adjacent  buses  may  be  linked  by  coupling  each  bus 
(one  symmetrically  and  one  antisymmetrically  to  prevent  direct  coupling  )  to  a  pair  of 
physical  ‘transfer’  qubits.  A  two  qubit  CN  gate  across  two  or  more  linked  buses  is 
accomplished  by  using  one  (two  logical-qubit)  Molmer-Sorensen  gate  per  bus  transfer, 
along  with  single  qubit  gates  for  the  control  and  target  qubits.  Linking  buses  in  this 
mannner  thus  allows  for  the  construction  of  a  quantum  computer  with  the  topology  of  a 
scalable  open  network. 
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Quite  recently  it  has  been  evidenced  that  quantum  states  may  support  new  kind  of  orders, 
which  cannot  be  characterized  by  broken  symmetries  and,  thus,  cannot  be  described  by  the 
conventional  Ginzburg-Landau  theory  [1],  Quantum  orders  may  be  viewed  as  the  pertinent 
description  of  the  pattern  of  the  quantum  entanglement  in  a  quantum  many-body  ground-state. 
Quantum  or  topological  orders  have  already  been  extensively  used  in  the  analysis  of 
Fractional  Quantum  Hall  systems  [2],  leading  to  an  elegant  explanation  of  their  robust 
(against  a  weak  but,  otherwise,  arbitrary  perturbation)  ground-state  degeneracy  and 
evidencing  the  intimate  connection  between  the  ground-state  and  the  statistics  of  quasi¬ 
particles.  The  robustness  of  the  ground-state  degeneracy  makes  quantum  ordered  states  very 
relevant  candidates  for  the  engineering  of  quantum  devices  naturally  taming  the  intrinsic 
decoherence  of  other  solid  state  devices  [3],  Very  recent  studies  hint  to  the  new  and  exciting 
possibility  that  the  topology  of  Josephson  junction  networks  (JJN)  may  be  crucial  for 
inducing  novel  and  unexpected  macroscopic  quantum  phenomena  [4]  as  well  as  for  opening 
[5]  to  the  possibility  for  an  explicit  realization  of  a  quantum  ordered  state. 

In  this  talk  I  shall  focus  the  attention  on  the  recently  discovered  possibility  that  JJN  built  on 
suitable  graphs  may  support  new  and  interesting  quantum  macroscopic  states  induced  only  by 
a  pertinent  engineering  of  the  geometry  and  topology  of  the  graph  supporting  the  network  [4]. 
In  particular,  I  shall  consider  JJN  built  on  comb-like  graphs  and  star-graphs  evidencing  the 
fact  that  -  already  for  this  very  simple  graph  topology  -  the  array  exhibits  quantum 
macroscopic  coherence  at  low  temperatures.  After  discussing  the  theoretical  possibility  of  a 
spatial  Bose-Einstein  condensation  (BEC)  of  Cooper  pairs  for  a  JJN  built  on  a  comb-graph,  I 
also  shall  discuss  here  a  simple  experiment  which  could  enable  to  observe  BEC  in  these 
systems  . 

A  network  of  superconducting  Josephson  junctions  is  usually  described  by 

H  =  -EJ'Zcos^l-^l)+EcY,nl 

(w)  '■ 

On  each  site  i  of  the  network  there  is  a  superconducting  grain  and  the  junctions  are  located 
between  neighboring  sites  with  Josephson  energy  Ej ;  Ec-  2e2/C0 ,  where  C0  is  the  grain 
self-capacitance  and  e  the  electron  charge,  is  the  charging  energy.  $  is  the  phase  of  the 
superconducting  order  parameter  at  site  labeled  by  i  and  its  conjugate  variable  n 
([$,«,]  =  <5y)  describes  the  number  of  Cooper  pairs  in  the  i-th  superconducting  grain.  The 
parameter  E,  is  related  to  the  critical  Josephson  current  7(0)  by  the  relation  7(0)  -leEjtti. 
The  summation  runs  over  all  the  distinct  nearest-neighbors  pairs.  By  properly 

engineering  the  insulating  support  on  which  the  superconducting  grain  are  placed,  one  can 
realize  star-shaped  or  comb-shaped  networks:  in  Fig.  1  a  schematic  description  of  a  Josephson 
junction  comb  network  is  reported. 

The  computation  of  the  Cooper  pairs  number  for  comb-  [6]  and  star-network  [7]  yields  that, 
below  the  critical  temperature  Tc  at  which  BEC  occurs,  the  critical  Josephson  current  IC(T) 
of  a  junction  far  away  from  the  backbone  (center)  of  the  comb  (star)  is  given,  in  the  limit 
Ec  ->0,  by 


r(°)  t 

1  *  C 

where  7(0)  is  the  critical  current  of  the  same  Josephson  junction  for  T>TC  (see  inset  of 

Fig.l).  This  suggests  that  the  measurement  of  the  critical  current  of  a  suitably  chosen  junction 
should  enable  to  detect  the  effect  of  the  topology  induced  BEC  in  graph-shaped  network. 
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Fig.l.  A  schematic  plot  of  a  typical  experimental  setup  with  50  x  50  Josephson  junctions  placed  at  the  sites  of 
the  comb-shaped  graph.  Notice  that  superconducting  grains  belonging  to  neighbouring  vertical  fingers  of  the 

comb  are  not  weakly  coupled.  Inset:  Predicted  Josephson  critical  current  (in  units  of  7(0))  vs.  TITC  of  a 
junction  located  far  away  from  the  horizontal  backbone. 
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Very  recently,  Gorelik  et  al.  [1]  proposed  a  very  appealing  setup,  the  Cooper  pair  shuttle,  able  to  create  and 
maintain  phase  coherence  between  two  distant  superconductors.  It  consists  in  a  superconducting  SET  transistor 
where  the  metallic  grain  is  driven  in  a  periodic  motion  between  the  two  electrods.  The  shuttling  effect  manifests  in 
the  fact  that  there  is  charge  transport  even  though  the  grain,  during  is  motion,  is  always  in  contact  with  only  one 
of  the  elecrodes.  Because  of  the  superconducting  electrodes,  the  shuttle  does  not  only  carry  charge,  as  in  the  normal 
metal  case,  but  it  also  establishes  phase  coherence  between  the  superconductors.  This  is  witnessed  by  the  presence  of 
a  steady  stste  Josephson  current. 

We  analyze  how  the  presence  of  the  environment  affects  the  coherent  transport,  both  in  the  average  current  and  its 
moments,  in  the  Cooper  pair  shuttle  [2],  For  the  environment  we  consider  the  case  when  this  is  determined  gate  voltage 
fluctuations.  The  expression  of  the  Josephson  current  depends  on  the  phases  of  the  superconducting  electrodes,  on 
the  dynamical  phases  acquired  during  the  motion  of  the  central  island  and  on  the  dephasing  rates.  Depending  on  the 
value  of  the  dynamical  and  superconducting  phases  the  critical  current  can  be  negative,  i.e.  the  system  can  behave  as 
a  7r-junction.  An  analysis  of  the  critical  current  as  a  function  of  the  dephasing  rates  reveals  another  interesting  aspect: 
the  Josephson  current  is  a  non-monotonous  function  of  the  damping  the  Josephson  current  can  increase  on  increasing 
the  coupling  to  the  environment.  The  typical  behaviour  of  the  Josephson  current  as  a  function  of  the  phases  is  shown 
in  Fig.l 

Cooper  pair  shuttling  is  a  result  of  a  non  equilibrium  steady  state  process.  Therefore,  to  better  characterize  the 
transport,  we  analyzed  supercurrent  noise  fluctuations.  This  should  be  contrasted  with  the  standard  Josephson  effect, 
where  the  supercurrent  is  an  equilibrium  property  of  the  system.  Finally,  we  briefly  discuss  the  finite  frequency 
spectrum  in  the  case  of  strong  dephasing.  Superimposed  to  the  peak  at  the  Josephson  energy,  there  are  oscillations  of 
frequency  of  the  order  of  period.  The  presence  of  these  oscillations  is  related  to  the  periodicity  of  the  island  motion. 
The  modification  of  these  fringes  as  a  function  of  the  phases  is  a  signature  of  the  coherence  in  the  Cooper  pair  shuttle. 

The  full-counting  statistics  was  analyzed  by  one  of  us  in  [3]. 

A  possible  experimental  test  of  our  results  which  does  not  require  any  mechanically  moving  part.  The  time  depen¬ 
dence  of  the  Josephson  couplings  can  be  regulated  by  a  time  dependent  magnetic  field  and  gate  voltage,  respectively. 
The  setup  consists  of  a  superconducting  nanocircuit  in  a  uniform  magnetic  field.  By  substituting  the  Josephson 
junction  by  SQUID  loops,  it  is  possible  to  control  the  Ej  by  tuning  the  applied  magnetic  field  piercing  the  loop.  A 
configuration  with  three  type  of  loops  with  different  area,  allows  to  achieve  the  require  time  dependence  by  means 
of  a  uniform  magnetic  field.  This  implementation  has  several  advantages.  It  allows  to  control  the  coupling  with  the 
environment  by  simply  varying  the  time  dependence  of  the  applied  magnetic  field.  The  time  scale  for  the  variation  of 
the  magnetic  field  should  be  controlled  at  the  same  level  as  it  is  done  in  the  implementation  of  Josephson  nanocircuits 
for  quantum  computation. 


FIG.  1:  Supercurrent  (in  units  of  e/T)  as  a  function  of  the  superconductor  phase  difference  <j)  and  of  the  phase  accumulated 
during  the  contact  to  one  of  the  electrodes  9. 


In  the  limit  in  which  the  Josephson  energy  is  much  larger  than  the  charging  energy,  the  Cooper  pair  shuttle 
implements  [4]  a  paradigm  model  for  quantum  chaos,  the  quantum  kicked  rotator.  In  this  limit  one  can  study 
dynamical  localization  and  by  further  coupling  the  shuttle  to  a  Cooper  pair  box  it  is  possible  to  measure  the  dynamical 
fidelity  as  suggested  in  [5]. 
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Controlling  the  dynamics  of  a  complex  quantum  system  is  at  the  heart  of  quantum  information  processing  [1], 
However  in  any  real  device  the  computational  degrees  of  freedom  entangle  with  the  environment  leading  to  decoher¬ 
ence  [2].  Open-loop  Bang-Bang  (BB)  quantum  control  techniques  have  been  proposed  as  a  way  to  achieve  an  effective 
decoupling  from  the  environment  [3,  4].  They  may  be  operated  by  a  sequence  of  very  strong  external  pulses  separated 
by  a  time  At  [3].  In  the  limit  At  — >  0  full  decoupling  [3,  4]  of  the  unwanted  interactions  is  achieved  since  the  effective 
dynamics  produced  by  them  is  averaged  out  by  BB.  The  physics  in  this  limit  is  a  manifestation  of  the  quantum  Zeno 
effect  [3,  5].  In  practice  At  is  finite  especially  if  full-power  pulses  are  generated.  Thus  decoupling  is  not  perfect  but 
the  Zeno  limit  physics  is  still  a  good  description  if  At  <C  7_1,  the  typical  time  scale  of  the  environment  [3,  4].  If  7  is 
large  one  may  argue  that  BB  chops  noise  and  frequencies  w  <  1/At  <  7  are  averaged  out.  This  optimistic  scenario 
could  foresee  important  applications  to  solid  state  coherent  devices,  where  low-frequency  noise  [6]  is  one  of  the  major 
problems  for  quantum  state  processing  [7-9]  and  has  been  exploited  in  a  series  of  recent  papers  [10-12], 

Here  we  study  environments  of  dissipative  quantum  bistable  fluctuators  [7],  describing  charge  noise  due  impurities 
close  to  a  solid-state  qubit  [7,  8,  13]  which  are  responsible  of  a  distictive  phenomenology  recently  observed  [14].  We 
consider  a  qubit  (' Hq  =  -§  oz  -  ^  crx)  coupled  to  an  impurity.  The  Hamiltonian  has  the  general  form  [12] 

H  =  Hq  -  i  E  +  He  +  V(f)  (1) 

where  He  describes  an  impurity  level  occupied  by  a  localized  electron,  tunneling  to  a  fermionic  band  whith  switching 
rate  7.  The  charge  in  the  impurity  is  coupled  to  the  qubit,  E  —  vnimp,  and  operates  as  a  stray  random  bias  of 
amplitude  v.  Control  is  operated  by  V(f),  representing  a  sequence  of  instantaneuos  7r-pulses  about  the  x  axis  of  the 
Bloch  sphere.  The  evolution  operator  of  the  Hamiltonian  (1)  is  [<rx<S]27V,  where  S  =  exp(-iHAt)  with  V(t)  =  0  is  the 
evolution  between  pulses.  In  order  to  take  into  account  memory  effects  of  the  environment  we  denote  with  p(t)  the 
reduced  density  matrix  (RDM)  of  the  qubit  plus  localized  level  and  use  a  standard  master  equation  to  find  the  map 
£ai[p(i)]  =  p{t  +  At)  =  eCAtp(t).  After  t  =  2NAI  we  get  p(t)  as  [P  ■  eCAt]2N  p(Q)  where  V  is  the  superoperator  of  the 
pulses.  The  RDM  of  the  qubit  pQ  =  Tr<j[p(t)]  is  obtained  by  tracing  out  the  localized  level  at  the  end  of  the  protocol. 


FIG.  1:  Left  panel  -  For  fixed  t  —  IO7-1  we  plot  T N{t)/g2  for  BB  procedures  with  N  echoes.  The  parameter  is  p  =  v/'y. 
A  gaussian  environment  with  the  same  power  spectrum  would  give,  for  arbitrary  g,  the  curve  here  labeled  with  g  =  0.1.  In 
the  inset  F(t)  for  g  —  1.1  without  (thick  line)  and  with  pulses  compared  with  results  of  the  stochastic  Schrodinger  equation 
(crosses). 

Central  panel  -  We  plot  S  =  ln[Tr(p®)2]  at  fixed  t  =  87“ 1  for  protocols  with  N  echo  pair  of  pulses.  Curves  are  parametrized 
by  g  =  (fli  —  f2o)/7.  Here  e  =  A,  v/Q.  =  0.2.  In  the  inset  it  is  shown  that  S  scales  as  g2  for  /V  >  1. 

Right  panel  -  BB  control  of  1//  noise  for  A  =  0.  The  analytic  r*r(t)  at  t  =  2N  At  (symbols)  is  compared  with  the  evolution 
with  no  pulses  (thick  solid  line).  Noise  is  generated  with  Np  fluctuators  with  rates  7;  distributed  from  10 *  Hz  to  10 10  Hz. 
Noise  level  oc  v2 N ji  is  fixed  at  a  value  typical  of  experiments  in  charge  qubits:  it  is  realized  with  coupling  v  =  9.23  •  109  Hz,  for 
9  =  6  (full  triangles),  7  (circles)  and  8  (squares),  with  Np  =  6  •  1017_2?.  Crosses  are  stochastic  Schrodinger  simulations  with 
10s  realizations  of  the  1//  process,  for  9  =  7. 


The  major  step  of  this  program  is  the  diagonalization  of  the  8x8  supermatrix  [V  ■  e£At]  . 

In  some  cases  the  calculation  can  be  carried  out  analytically,  as  for  A  =  0  [12].  We  find  that  decoupling  of  this 
environment  is  sensitive  to  details  of  its  dynamics.  If  pulses  are  not  very  frequent  it  shows  a  rich  variety  of  behaviors, 
suggesting  that  BB  may  also  be  used  for  spectroscopy  of  a  solid-state  environment.  An  important  parameter  is 
g  =  (0X  —  flo)/7,  where  ft*  are  the  splittings  of  the  qubit  when  the  fluctuator  sits  in  state  i,  which  quantifies  the 
non-gaussianity  of  the  environment  [7].  Slow  flnctuators  (g  >  1)  are  responsible  for  inhomogeneuos  broadening  and 
for  a  physics  dominated  by  memory  effects  where  decoherence  depends  strongly  on  details  of  the  protocols.  We  study 
the  decay  of  the  coherences  of  the  qubit,  given  by  T^(t)  =  In  J  p°bb(0)+pIa(ol  (Fig-  1  left).  At  any  fixed  t  =  27V At, 
|F;v(t)|  decreases  monotonically  when  the  pulse  frequency  1/At  increases,  which  shows  that  BB  effectively  suppresses 
RTN  and  in  this  regime  |rV(t)|  ~  g2  shows  universal  behavior.  On  other  hand  for  27V  <  7 1  qualitative  differences  in 
the  behavior  are  apparent  for  g  <  1  and  g  >  1.  The  physics  for  A  ^  0  is  even  richer.  We  study  the  degree  of  purity 
S  =  lnTr(p<3)2,  efficient  decoupling  corresponding  to  S  =  0,  i.e.  localization  in  a  “Zeno  subspace”  [5]  which  is  a  pure 
state.  Results  (Fig.  1  central)  show  that  for  frequent  pulses  S  «  0,  but  decoupling  of  a  slow  fluctuator,  g  >  1,  requires 
a  very  large  TV.  Again  we  find  universal  behavior,  S  ~  g2  (inset  of  Fig.  1  central).  Instead  for  a  smaller  number  of 
pulses  it  may  happen,  especially  for  g  >  1  that  S  is  not  monotonic  with  TV,  including  the  possibility  that  the  qubit 
decays  faster  than  in  absence  of  pulses  [3],  reminiscent  of  the  inverse  Zeno  effect  [5]. 

In  order  to  treat  1//  noise  we  now  extend  our  formalism  to  a  “multi-mode”  environment.  We  first  generalize  the 
results  at  A  =  0  of  Ref.  [3],  to  an  arbitrary  (non-gaussian)  environment.  For  the  evolution  between  two  pulses  at  t\  and 
*2  we  can  use  [H,  az]  =  0  and  following  the  same  steps  of  Ref.  [3]  we  can  write  the  evolution  operator  for  a  BB  procedure 
at  t  =  27V At  as  S'2N(t)  —  [e'I 11*KE+5',‘^,Ate“'(?<«“5‘'»®)A,j jV.  In  the  overall  BB  procedure  az  is  conserved,  so  we 
need  only  off  diagonal  entries  of  the  RDM  of  the  qubit,  in  the  az  basis  p®ai(t)  =  p®a,  (0)  Trg  (t\a)  u^vS^Ict')} 

where  we  assumed  factorized  initial  conditions.  Here  <S2iv(i|cr)  =  (cr  |tS2Jv  (t)  ]  rr)  expresses  the  conditional  evolution 
of  the  environment  under  a  well  defined  sequence  of  a  —  ±1.  The  trace  factorizes  if  the  environment  is  made  of 
noninteracting  “modes”,  so  we  can  apply  this  result  to  a  set  of  noninteracting  impurities  originating  1//  noise.  In 

|  _(*>)  1  _(»?)  /£\  j 

particular  the  decay  of  the  coherences  is  V pj(t)  =  In  pi{o)+pd'(g)  ■  Results  in  the  right  panel  of  Fig.  1  show  that 
frequent  pulses  (curves  with  many  symbols)  drastically  reduce  decoherence.  For  noise  levels  typical  of  experiments 
with  charge  qubits  (Fig.  1  right)  the  pulse  rate  for  substantial  recovery  is  practically  independent  on  v.  The  situation 
may  change  if  a  broad  distribution  of  the  couplings  v  is  considered  [10].  The  physics  is  richer  for  A  ^  0  and  possibly 
compensation  of  1//  noise  is  non  monotonic  for  increasing  pulse  frequency,  as  it  happens  for  a  single  impurity. 
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Since  manipulations  of  qubits  are  constrained  by  the  quantum  mechanical  rules, 
several  classical  information  tasks  can  not  be  perfectly  extended  to  the  quantum  world. 
The  more  relevant  limitations  in  quantum  information  processing  are  the  impossibility 
to  perfectly  clone  any  unknown  qubit  |  <|)  >  and  to  map  it  in  its  orthogonal  state  |  (j)1  >. 
Even  if  these  two  processes  are  unrealizable  in  their  exact  forms,  they  can  be  optimally 
approximated  by  the  so-called  universal  quantum  machines:  the  optimal  quantum 
cloning  machine  (UOQCM)  and  the  universal-NOT  (U-NOT)  gate,  which  exhibit  the 
minimum  possible  noise.  Investigation  of  these  optimal  transformations  is  important 
since  it  reveals  bounds  on  optimal  manipulations  of  information  with  quantum  systems. 

The  U-NOT  gate  and  the  UOQCM,  have  been  contextually  realized  by  adopting 
the  process  of  stimulated  emission  generated  by  a  single  photon  into  an  optical 
parametric  amplifier  [1].  We  shall  present  how  it  is  possible  to  implement  the  1  -»  2 
UOQCM  and  the  1  ->  1  U-NOT  gate  by  slightly  modifying  the  quantum  state 
teleportation  protocol  [2].  In  this  case  the  UNOT  gate  is  transferred  in  a  different 
location,  we  deal  hence  with  the  teleportation  of  a  quantum  operation:  the  Tele-UNOT 
gate.  A  complete  experimental  realization  of  the  protocol  by  a  fully  linear  method  will 
be  discussed. 
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A  Monte  Carlo  method  which  allows  one  to  follow  each  electron  tunneling  event  has 
been  successfully  applied  to  simulate  transport  properties  of  a  single-electron  transistor 
(SET)  or  more  complicated  single  electronics  circuits.  This  method  gives  physical 
insight  into  the  processes  taking  place  in  the  simulated  system,  and  macroscopic 
ensemble  properties  can  be  calculated  using  large  ensembles  of  single  electron 
tunneling  events.  But  to  our  knowledge,  it  has  not  yet  been  formally  applied  to  quantum 
measurement  problems  by  a  SET  detector.  Here,  we  provide  such  an  investigation  and 
derive  the  quantum-jump  [1,2]  stochastic  master  equation  (or  quantum  trajectory 
equation)  for  a  single  superconducting  Cooper-pair-box  (SCB)  charge  qubit 
(generalization  to  other  charge  qubit  case  is  simple)  continuously  measured  by  a  SET 
[3],  We  also  show  that  the  master  equation  for  the  “partially”  reduced  density  matrix 
presented  in  Ref.  [3]  (referred  to  here  as  a  “partial”  course-grain  description)  can  be 
obtained  by  taking  a  “partial”  average  on  the  stochastic  master  equation  over  the  fine 
grained  measurement  records  of  the  tunneling  events  in  the  SET. 

If  a  measurement  is  made  on  the  qubit  system  and  the  results  are  available,  the  state  or 
density  matrix  is  conditioned  on  the  measurement  results.  If  the  subsequent  system 
evolution  after  the  measurement  is  concerned,  the  conditional,  stochastic  master 
equations  derived  here  should  be  employed.  A  set  of  typical  quantum  trajectories 
(stochastic  state  evolutions)  and  its  corresponding  measurement  record  are  shown  in  Fig. 
1  to  illustrate  the  conditional  dynamics.  We  consider  the  case  that  the  leading  tunneling 
process  in  the  SET  are  sequential  transitions  between  two  adjacent  charge  states  N=  0 
and  N+ 1=1  on  the  SET  island.  dNL/Rc  represents  the  number  (either  zero  or  one)  of 

electron  tunneling  events  seen  in  the  infinitesimal  time  dt  through  the  left  and  right 
junctions  of  the  SET,  respectively.  The  subscript  c  indicates  that  the  quantity  to  which  it 
is  attached  is  conditioned  on  previous  measurement  results.  As  expected,  due  to 
Coulomb  blockade,  the  SET  measurement  record  of  the  randomly  distributed  moments 
of  detections  shown  in  Fig.  l(i)  and  (j)  is  in  the  order  of  an  exactly  alternating  dNu  =  1 

and  dNRc  - 1  time  sequence.  The  joint  reduced  density  matrix  operator  of  the  SCB 
qubit  and  extra  charge  on  the  SET  island  can  be  written  as 
Pc  =  Pn,c  I  N)(,N  |  ~^Pn+-[  c  |  N  +  1>(A  -J- 1 1,  where  pNIN+l  c(t)  is  each  a  2x2  operator  in 
the  qubit  basis.  The  conditional  qubit  density  matrix  operator  alone  can  be  obtained  by 
writing  pc(t)  =  TrN[Rc(t)]  -  pNc(t)+ pN+lc(t) .  Figure  l(a)-(d)  show  the  time  evolution 

of  the  joint  SCB-qubit  and  SET-island  diagonal  density  matrix  elements  conditioned  on 
the  measurement  results  of  Fig.  l(i)  and  (f).  The  conditional  evolutions  of  the  qubit 
alone  shown  in  Fig.  1(e)  and  (I)  can  be  obtained  from  the  sum  of  the  joint  state 
evolutions  of  (a)  and  (b),  or  (c)  and  (d),  respectively.  The  probabilities, 

Pon,c  =  Trqb[Po/i,J  =  pZ,c  +  Pm,c  of  the  SET  island  state  alone  in  (g)  and  (h)  can  be 


obtained  by  summing  the  evolutions  in  (a)  and  (c),  or  (b)  and  (d),  respectively.  Each 
quantum  trajectory  mimics  a  single  history  of  the  qubit  state  in  a  single  run  of  the 
continuous  measurement  experiment.  The  conditional  evolutions  in  (a)-(h)  differ 
considerably  from  their  ensemble  average  counterparts. 

If  only  one  measurement  value  is  recorded  in  each  run  of  experiments  (e.g.,  the  number 
of  electrons  NR  that  have  tunneled  into  the  right  lead  or  the  drain  in  time  t)  and 

ensemble  average  properties  (e.g,  P(NR,t)  the  probability  distribution  of  finding  NR 

in  time  t )  are  studied  over  many  repeated  experiments,  the  quantum  trajectory  approach 
presented  here  will  give  the  same  result  as  the  master  equation  approach  of  the 
"partially"  reduced  density  matrix  [3],  This  is  demonstrated  in  Fig.  2.  However,  the 
possible  individual  realizations  of  quantum  trajectories  and  their  corresponding 
measurement  records  (e.g,  in  Fig.  1)  do  provide  insight  into,  and  aid  in  the  interpretation 
of  the  average  properties. 
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Figure  1.  Conditional  time  evolution  of  the  SCB  charge  qubit  measured  by  the  SET. 


Figure  2.  Probability  distribution 
P(NR,t  -  2000) .  Our  simulation  of  the 

normalized  histogram,  using  2000 
quantum  trajectories  and  their 
corresponding  detection  records,  is 
already  in  good  agreement  with  the  plot  in 
solid  line,  obtained  from  the  Fourier 
analysis  of  the  partially  reduced  density 
matrix  method  described  in  Ref.  [3]. 
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We  present  an  extensive  study  of  entanglement  properties  in  easy-plane  quantum  spin 
systems  with  a  quantum  phase  transition  driven  by  a  magnetic  field  applied  in  the  easy 
plane.  Making  use  of  quantum  Monte  Carlo  simulations,  we  are  able  to  monitor  the 
behavior  of  englement  in  the  ground  state  and  at  finite  temperature  as  a  function  of  the 
applied  field  for  different  lattice  geometries  of  the  system.  Our  calculations  focus  on  the 
entanglement  of  formation,  quantified  by  the  one-tangle  [1]  for  the  global  entanglement  of 
the  system  and  by  the  concurrence  [2]  for  the  pairwise  englement.  For  the  case  of  a  one¬ 
dimensional  spin  chain,  we  observe  that  the  entanglement  estimators  are  able  to  single  out 
with  high  precision  both  the  quantum  critical  point,  through  a  minimum  of  the  paiwise-to- 
global  entanglement  ratio,  and  the  known  occurrence  of  a  factorized  state  below  the 
critical  field,  through  the  vanishing  of  all  entanglement  estimators  [3],  The  vanishing  of 
entanglement  at  a  given  field  also  persists  at  finite-temperature,  so  that  the  system  displays 
a  temperature-resistant  entanglement  switch  effect.  We  then  extend  the  same  analysis  to 
the  case  of  a  spin  ladder  and  of  a  square  lattice,  and,  thanks  to  the  entanglement 
estimators,  we  find  that  the  existence  of  a  factorized  state  close  to  the  quantum  critical 
point  is  a  general  feature  independent  of  the  geometry. 
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Introduction 

The  reported  phase  qubit  is  dc  superconducting  quantum 
interferometer  (dc  SQUID)  containing  two  Josephson 
junctions,  which  are  characterized  by  a  second  harmonic  in 
their  current-phase  relations  (CPRs): 

lj  =  ICj  sin(  <p j)  +  7C/  sin(2  (pj),  j=  1,2,  ( 1 ) 

where  7C/  and  Icf  -  amplitudes  of  the  first  and  second 
harmonics  respectively,  (pj  -  Josephson  junction  phase.  This 
qubit  was  suggested  first  by  Charlebois  et  al.  [1]  to  be  realized 
on  the  base  of  high-Tc  superconductors.  The  advantages  of  this 
phase  qubit,  which  is  a  “quiet”  qubit,  are  the  absence  of  a 
state-dependent  spontaneous  current  in  the  loop  and  the 
uselessness  of  any  bias  magnetic  flux  [2], 

Qubit  Design  and  Properties 

The  total  Josephson  energy  of  the  interferometer  is  a  function  of  three  asymmetry 
parameters:  the  ratio  of  the  first  harmonic  amplitudes  of  the  junctions  r]  =  1C2 /  7c/  and 
the  algebraic  ratio  of  the  second-to-first  harmonic  amplitudes  in  the  current-phase 
relation  for  either  of  the  two  junctions  a\p  =  lc\p  / Ic\,i-  Normalized  inductance  of  the 
qubit  l  =  2ti7-7o/Oo  (O0  =  h/2e-  magnetic  flux  quantum)  must  be  as  small  as  possible  to 
eliminate  any  crosstalk  between  qubits.  When  one  considers  a  small  inductance  of  the 
qubit  (/  <<  1),  the  total  flux  O  through  the  loop  approaches  the  applied  magnetic  flux 
T>e.  In  this  case,  the  Josephson  energy  of  the  system  at  7c/  =  70  is  given  by 
E  1  1 

—  =  — (cos(—  +  O')  + — «,  cos {tp  +  26)  +  p  cos (—  -  6)  +  —  rja2  cos  fp  -  26)) ,  (2) 

El  2  2  2  2 

where  <p=(p  \-(p2=(p  e=2n  =  ( <p  i+  <p  2)  /2;  E\=  d>0  l0/2n-  characteristic  energy. 

When  no  external  flux  is  applied  (see  the  energy  profile  at  (p=  0  in  Fig.  I),  the 
energy  potential  is  always  symmetric  with  respect  to  6  =  (<p\+<p2)/2,  even  at  strong 
asymmetry  of  the  interferometer  ( |  a.\-a2  \  ~  1  or  |  rj-\  \  ~  1).  The  left  and  right  energy 
states  are  degenerate  and  form  a  basis  of  the  qubit.  This  qubit  is  “quiet”  because  no 
external  field  and  no  bias  current  are  required  to  form  the  qubit  states.  At  the  same  time 
there  is  threshold  value  a*  =  -0.6  of  the  normalized  amplitudes  of  second  harmonic, 
that  is  critical  for  the  phase  qubit  existence.  If  either  oq  or  a2  is  larger  then  a*,  the 
interferometer  energy  does  not  show  any  double-well  potential  profile. 

If  external  magnetic  field  is  applied,  the  potential  wells  are  lifted  and  the  symmetry 
of  the  potential  becomes  broken  at  rj  f  1  and  a\  f  a2  (unequal  junctions).  But  in  the 
extreme  case  of  inequality  of  the  junctions  (p  «  1  or  tj  »  1)  meaning  transformation 
of  the  two-junction  interferometer  to  a  one-junction  interferometer,  the  energy  potential 
is  symmetric  again  and  the  basis  states  are  degenerate  at  any  magnetic  fields  (at  least  up 
to  <pe~  2n). 


Fig.  1.  Josephson  energy  of  the  qubit  (cq  =  a2  =  -1,  rj  =  0.8)  versus  ©  and  cp  =  <pe  (left)  as  well  as  one- 
dimensional  profile  (right)  of  the  energy  for  different  values  of  the  normalized  magnetic  flux  <pe  applied. 

Evaluation  of  the  Second  Harmonic  Amplitude 

There  is  an  experimental  evidence  that  the  second  harmonic  can  reveal  itself  by 
formation  of  an  additional  side  peak  in  the  dependence  of  the  dc  SQUID  critical  current 
on  applied  magnetic  field  Ic(H)  [3]-[4],  The  second  harmonic  impact  on  the  dependence 
7c(®e)  has  been  studied  in  detail  numerically  for  different  values  of  the  normalized 
interferometer  inductance  /  by  means  of  computational  solution  of  the  simultaneous 
equations  (1)  and  equation  for  the  junction  phase  relation 

cp\-  <pi=  (p&~  ( l\i\  -  hh)  (3) 

to  find  a  maximum  value  of  the  superconducting  transport  current  through  the 
interferometer  ie  =  i\+i2  as  a  function  of  the  applied  magnetic  flux  <pc  (here  ie  -  1J  Io, 
i\,2  -  I\,2  /  Io,  l\,2  -  normalized  arm  inductances). 

It  has  been  found  that  at  Ic\  ~  Io  and  Ici*  =  0  the  additional  side  peak  caused  by 
Ic2JI  decreases  rapidly  with  inductance  /  of  the  interferometer  and  vanishes  completely  at 
/»  0.8  [5],  If  />  1,  this  peak  can  exist  only  at  sufficiently  small  values  of  the  first 
harmonic  amplitude  Ici .  The  dc  SQUID  inductance  /crit,  which  is  critical  for  the  side 
peak  existence  in  the  dependence  /c(<&e),  was  calculated  numerically  as  a  function  of 
both  the  first  and  second  harmonic  amplitudes. 

The  obtained  results  for  the  SQUID  critical  current  dependence  on  applied 
magnetic  field  allow  evaluation  of  the  second  harmonic  amplitude  from  the  side  peak 
parameters.  This  is  very  important  for  realization  and  optimization  of  the  reported 
“quiet”  phase  qubit. 
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We  present  a  novel  scheme  for  performing  a  conditional  phase  (CPHASE)  gate 
between  two  spin  qubits  in  adjacent  semiconductor  quantum  dots  through  delocalized 
single  exciton  states  [1].  The  delocalized  states  sire  formed  through  the  resonant 
energy  transfer  or  Foerster  mechanism,  in  which  an  exciton  can  transfer  between  two 
resonant  dots  via  an  intermediate  virtual  photon. 

We  therefore  consider  two  resonant  quantum  dots,  each  containing  a  single  excess 
conduction  band  electron  whose  spin  embodies  the  qubit.  Excitons  may  be  selectively 
created  depending  on  the  spin  state  of  the  dot  through  the  Pauli  blocking  effect,  in 
which  a  single  circularly  polarized  laser  pulse  is  applied,  resonant  with  the  s-shell 
heavy-hole  exciton  creation  energies.  These  excitons  are  necessarily  created  in  the 
z- angular  momentum  states  given  by  \3/2hh,  — l/2e).  and,  as  noted  in  Refs.  [[3,  2]], 
such  a  polarized  laser  will  only  create  an  exciton  on  a  quantum  dot  if  its  excess 
electron  is  in  the  spin  state  \mz  =  1/2),  due  to  the  Pauli  exclusion  principle. 

Since  the  Foerster  interaction  is  non-magnetic,  the  single  exciton  created  from  the 
1 11)  state  is  delocalized  and  has  a  different  energy  from  the  localized  single  excitons 
created  from  the  |10)  and  (01)  states  (see  Fig.  1).  It  is  therefore  possible  to  selectively 
create  an  exciton  from  the  |11)  state  only,  and  by  applying  a  27r  pulse,  a  phase  of 
— 1  may  be  picked  up  on  this  state,  whilst  the  other  states  remain  unchanged  (see 
Fig.  2).  This  is  the  CPHASE  gate  which  is  a  universal  two  qubit  gate  for  quantum 
computing. 

Following  on  from  work  in  Refs.  [3,  4],  we  also  show  that  arbitrary  single-qubit 
rotations  may  be  realized  to  a  high  fidelity  with  current  semiconductor  and  laser  tech¬ 
nology.  Two  distinct  single  qubit  operations  are  required  for  this.  First,  a  Z-rotation 
gate  is  achieved  by  again  exciting  single  excitons  selctively  using  Pauli  blocking.  Sec¬ 
ond,  an  A -rotation  gate  is  achieved  by  performing  Raman  transitions  between  the 
two  qubit  states,  through  the  light  hole  exciton  (see  Fig.  3). 
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Figure  1:  Schematic  of  the  energy  level  structure  of  two  resonantly  coupled  quantum 
dots,  showing  its  dependence  on  the  spin  state  of  the  excess  electron  on  each  dot. 


Figure  2:  Amplitude  (left)  and  phase  (right)  of  states  |10)  and  |11)  during  the 
CPHASE  gate  operation,  with  Foerster  stragnth  VF  =  0.85  meV,  biexciton  bind¬ 
ing  energy  VXx  =  5  meV,  and  exciton  creation  energy  tja  =  2  eV.  Curves  are  shown 
for  a  variety  of  Rabi  frequencies  fi. 


Figure  3:  Population  of  the  state  of  a  single  qubit  during  the  proposed  A-rotation 
gate  induced  by  a  two-laser  Raman  transition.  Various  values  of  the  exciton  decay 
rate  7  and  laser  detuning  v  are  shown  with  uia  =  2  eV,  and  fl  =  1.33  meV  for  both 
lasers. 
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Recently  amazing  progress  has  been  made  in  the  achievement  of  increasingly 
degenerate  regimes  of  trapped  atomic  Fermi  gases.  After  reaching  the  degenerate 
regime,  experiments  have  dealt  with  mixtures  of  two  hyperfine  levels.  This  allows 
quite  conveniently  to  adjust  the  effective  interaction  between  atoms,  employing  a 
Feshbach  resonance.  At  low  temperatures  the  scattering  between  atoms  is 
characterized  by  their  s-wave  scattering  length.  By  slowly  varying  the  magnetic 
field  one  can  sweep  through  the  magnetic  Feshbach  resonance.  Starting  with  a  weak 
attractive  interaction  between  atoms  of  different  hyperfine  levels,  i.e.,  small 
negative  scattering  length,  making  it  more  and  more  negative,  and  finally  have  it 
jump  to  a  positive  value  at  the  other  side  of  the  resonance,  corresponding  to  a 
repulsive  interaction  between  atoms  leading  to  a  formation  of  diatomic  molecules. 
Utilizing  such  magnetically  tuned  scattering  resonances  between  two-component 
Fermi  gases  experimentalists  have  been  able  to  produce,  Bose-Einstein  condensates 
of  molecules  formed  by  Fermionic  atoms  of  6Li  and  40K  as  well  as  to  probe  the, 
theoretically  predicted  crossover  from  a  Bose-Einstein  condensate  (BEC)  to  a 
Bardeen-Cooper-Schriefer  (BCS)  superfluid . 

An  interesting  feature  of  a  couple  of  such  Feshbach  sweep  experiments  [1,2]  is  that 
the  experimental  efficiency  observed  in  the  conversion  of  ultracold  Fermi  atomic 
gases  of  40K  and  6Li  atoms  into  diatomic  Bose  gases  was  limited  to  0.5.  By 
constructing  [3]  the  many-body  state  we  show  that  the  initial  state  preparation 
performed  in  such  experiments  involving  two-component  Fermi  gases  yields  a 
mixture  of  even  and  odd  parity  pair-states,  composed  of  two  subsystems 
corresponding  to  two  different  spin  states.  Through  the  symmetry  of  the  initial 
many-body  state  we  are  able  to  explain  the  50%  limited  transfer  efficiency  in 
transforming  atoms  into  molecules  obtained  in  the  above  experiments  as  well  as  the 
fact  that  the  Feshbach  resonance  sweep  does  not  yield  a  swift  heating  of  the 
remaining  atomic  Fermionic  gas  as  was  theoretically  predicted  [4], 

In  very  low  temperatures  the  Feshbach  resonance  is  due  to  s-wave  interactions, 
therefore  the  scattering  atoms  need  to  have  a  spatially  symmetric  wave  function  in 
order  to  create  a  molecule.  In  considering  pairs  of  atoms  in  the  initial  state  of  the 
system,  one  from  each  spin  state,  half  the  atomic  pairs  are  anti-symmetric  spin- 
states  and  half  are  spin  symmetric.  More  explicitly,  the  reduced  two-particle  density 
matrix  obtained  by  tracing  out  all  but  one  particle  of  one  spin  state  and  another 
particle  of  the  other  spin  state  contains  50%  spin  anti- symmetric  spin-states, 
interacting  via  an  s-wave  (even  parity  spatial  state),  and  50%  spin  symmetric  spin 
superposition  which  are  effectively  non-interacting  at  such  low  temperature.  Thus 
when  sweeping  through  the  resonance  sufficiently  slow  so  that  the  Landau-Zener 
transition  is  traversed  adiabatically  it  is  clear  that,  since  only  even  parity  states  can 
produce  molecules  through  the  Feshbach  sweep,  the  conversion  of  ultracold  Fermi 
gases  into  a  diatomic  Bose  gas  is  limited  to  0.5. 


To  obtain  a  transfer  efficiency  which  is  greater  than  0.5,  as  is  commonly  achieved 
in  current  experiments,  the  odd  parity  spin-symmetric  states  must  decorrelate  before 
the  constituent  atoms  can  interact  again  via  the  Feshbach  resonance.  Therefore  in 
achieving  a  higher  efficiency  the  sweep  rate  was  much  slower,  approaching  "close 
to  thermal  equilibrium"  conditions.  Moreover  in  the  relatively  rapid  sweep 
experiments  [1,2]  in  which  the  Feshbach  resonance  sweep  rate  is  sufficiently  slow 
to  pass  adiabatically  through  the  Landau  Zener  transition  but  faster  than  "the 
collision  rate"  in  the  gas,one  expects  a  rapid  heating  of  the  Fermionic  gas.  However 
even  though  the  Feshbach  sweep  produces  holes  in  the  Fermionic  atomic  sea,  these 
holes  can  not  be  filled  by  the  remaining  atomic  population,  i.e.,  the  population  of 
atoms  which  have  not  been  converted  to  molecules,  since  their  composing  atoms 
are  non-interacting,  due  to  their  spatial  symmetry  which  controls  their  interaction  at 
ultra-low  temperatures.  This  explains  the  experimental  observation  that  the 
remaining  atoms  do  not  heat. 
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Linear-response  conductance  of  the 
normal  conducting  single-electron  pump 
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For  single-electron  devices  the  model  Hamiltonian  matches  its  experimental 
realization  very  closely.  This  has  been  proven  most  impressively  for  the  normal¬ 
conducting  single-electron  transistor  (SET)  with  the  aid  of  Quantum-Monte-Carlo 
methods  [1,2].  However,  the  lowest  order  perturbation  theory,  the  so  called  sequential 
model,  gives  only  a  qualitative  description  in  most  practical  cases.  It  ignores  quantum 
fluctuations  due  to  the  coupling  of  charge  states  by  the  finite  conductance  of  the  tunneling 
contacts,  which  are  fundamental  constituents  of  the  devices.  Perturbation  expansion  (PE) 
has  been  suggested  do  deal  with  the  impact  of  the  coupling.  Especially  if  the  Coulomb 
blockade  is  lifted  by  lining  up  the  energy  of  two  neighboring  charge  states  a  situation 
arises  which  calls  for  sophisticated  treatments.  PE  methods  developed  by  Scholler  and 
Konig  [3,4]  as  well  as  by  Grabert  and  Goppert  [1,5]  have  been  proven  to  describe  the 
SET  at  the  degeneracy  point  quantitatively  up  to  G  =  1.5  e2/h  [2,6],  Despite  this  success 
doubts  on  the  applicability  of  PE  in  more  general  cases  has  occurred  lately.  Even  for  the 
SET,  the  validity  range  away  from  degeneracy  is  not  clear.  E.g.  in  the  blockaded  regime 
deviation  between  PE  and  experimental  results  are  stronger  than  at  the  degeneracy  point 
while  the  quantum  Monte-Carlo  data  are  in  full  accordance  with  the  experiments. 
Furthermore,  intermediate  expressions  of  the  PE  calculation  are  burdened  by  cut-off 
parameters.  In  many  expressions  for  the  SET  -  including  the  mean  charge  on  the  transistor 
island  and  the  conductance  -  these  cut-offs  cancel.  But  this  happens  to  be  somewhat 
fortunate,  since  no  general  rule  is  known  which  guarantees  the  cancellation.  E.  g.  the 
occupation  probability  of  the  individual  charge  states  contain  the  cut-off  parameter 
explicitly. 

To  get  more  insight  into  the  applicability  of  PE  high  quality  conductance  data  on 
single-electron  devices  other  than  the  SET  are  desirable.  Here  we  present  measurements 
of  the  linear  response  conductance  of  single-electron  pumps  (SEP).  We  have  fabricated 
single-electron  pumps  with  aluminum  island  and  aluminum-oxide  tunnel  barriers  in 
shadow  evaporation  technique.  Usually  single-electron  pumps  consists  of  a  serial 
arrangement  of  three  tunneling  contacts  forming  two  island;  the  electrostatic  potential  on 
each  island  can  be  adjusted  via  gate  electrodes.  However,  in  this  simple  layout  only  the 
serial  conductance  of  the  contacts  can  be  measured  and  the  distribution  among  the 
individual  contacts  remains  unknown.  To  make  the  determination  of  the  conductance  of 
each  contact  possible,  we  replaced  each  outer  contacts  by  two  contacts.  In  the  resulting 
device  conductance  parameters  can  be  measured  along  different  current  paths.  From  a  set 
of  these  parameters  the  individual  conductances  can  be  determined.  In  the  final 
experiments  the  two  outer  contacts  on  the  left  and  the  right  island  are  connected  in 
parallel  to  the  source  and  drain  voltage  source  and  act  in  the  same  manner  as  a  single 
contact. 


We  report  the  measurements  on  three  well  characterized  samples  with  varying  asym¬ 
metry  of  the  conductance  between  inner  and  outer  tunneling  contacts.  In  all  cases  the 
conductance  is  of  the  order  of  the  conductance  quantum  e1 2 3 4 5 6/h.  This  results  in  strong 
quantum  fluctuations  of  the  charge  state,  which  manifest  themselves  -  similar  as  in  the 
case  of  the  SET  -  by  a  logarithmic  contribution  to  the  linear-response  conductance  at  low 
temperatures. 

In  addition  in  an  easily  accessible  parameter  regime  we  find  a  non-monotonic 
behavior  of  the  conductance  on  temperature,  which  is  correctly  described  by  the 
sequential  model.  Phenomena  of  this  kind  are  to  be  expected  in  all  single-electron  devices 
which  are  more  complex  than  the  SET.  They  can  uncover  internal  characteristics 
unaccessible  by  other  means  -  here  e.g.  the  asymmetry  of  the  conductance  between  the 
inner  and  outer  tunneling  contacts  which  is  not  directly  measurable  in  the  most  natural 
SEP  layout. 


Linear  response  conductance  of  three  SEP  as  a 
function  of  temperature.  Shown  are  the 
minimal  (o)  and  the  maximal  (•) 
conductance,  as  well  as  the  conductance  at  a 
highly  symmetric  point  of  the  SEP  gate 
voltages  (0).  The  dashed  lines  are  the  result  of 
an  analytical  solution  of  the  sequential  model 
obtained  by  restricting  the  model  to  four 
charge  states.  As  solid  lines  the  outcome  of  the 
full  sequential  model  with  independently 
determined  sample  parameters  is  shown.  For 
the  thin  solid  line  a  term  of  the  form  a 
log(kBT/Eco)  has  been  added  to  the  maximal 
conductance  as  calculated  from  the  sequential 
model.  Here  a  is  a  fitting  parameter  and  EOT 
measures  the  capacitive  coupling  between  the 
two  SEP  islands. 
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Entanglement  Assisted  Metrology 
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We  present  a  new  approach  to  the  measurement  of  a  single  spin  state,  based  on  NMR 
techniques  and  inspired  by  the  coherent  control  over  many-body  systems  envisaged  by 
Quantum  Information  Processing. 

We  propose  to  use  the  interaction  of  the  spin  of  interest  with  an  ensemble  of  spins  (that 
we  call  Spin  Amplifier)  in  order  to  enhance  the  signal  from  the  single  spin,  up  to  the 
point  where  it  is  detectable  by  the  usual  inductive  means.  This  method  demands  to 
control  the  polarization  of  about  106  spins,  by  inducing  a  coherent  dynamics  modulated 
by  the  interaction  with  a  single  spin. 

The  ensemble  of  spins  forming  the  Amplifier  is  first  prepared  in  a  highly  polarized 
state;  it  is  then  put  in  contact  with  the  target  spin  (the  spin  we  want  to  detect),  which  has 
already  collapsed  into  one  of  the  two  states  |0)  or  |1). 

A  simple  quantum  circuit,  which  illustrates  how  a  collective  measurement  can  provide 
knowledge  of  the  single  spin  state,  consists  of  a  train  of  Controlled  Not  (CNot)  gates 
between  the  target  spin  and  each  of  the  Amplifier  spins:  These  gates  effectively  copy 
the  state  of  the  single  spin  onto  the  state  of  the  Amplifier  spins.  At  the  end,  the 
measurement  of  the  total  polarization  in  the  Amplifier  (P  ocEj  (\|/|a‘z|\|/»  will  indicate  the 
state  of  the  target  spin.  This  first  scheme  demands  that  the  ensemble  spins  are 
independently  addressable,  and  the  interaction  with  the  single  spin  can  only  cause  local 
perturbations  on  single  spins. 

Instead  of  demanding  addressability  of  individual  spin,  we  can  develop  equivalent 
schemes  that  connect  better  to  the  control  we  have  on  the  physical  system  (and  hence 
realizable  in  the  near  term),  relying  only  on  the  collective  behavior  of  the  system  and  on 
entanglement  among  the  Amplifier  spins.  With  the  creation  of  a  macroscopic  entangled 
state  [1]  (cat-state)  in  the  Amplifier,  the  device  acts  as  a  collective  single  spin  and  thus 
is  globally  affected  by  a  sole  interaction  with  the  target  spin. 


COOH),  while  the  3  Carbons  compose  the  Amplifier.  The  spins  interact  via  the  scalar 
(J-)coupling  IijCTz'azJ. 

Before  applying  the  circuit  of  Fig.  1,  we  first  put  the  Proton  spin  into  the  identity 
state  and  prepare  the  Carbons  in  the  pseudo-pure  state  jOOO)  [2].  The  state  of  the  Proton 
is  a  superposition  of  the  two  decohered  states  |0)(0|  and  |1)(1|:  We  can  therefore 
effectively  perform  two  experiments  in  parallel,  and  read  out  the  outcomes  from  just 
one  spectrum,  where  the  lines  are  separated  by  the  Carbon-Proton  coupling  splitting. 


FIG.  2:  Experimental  Results: 

a)  Initial  State:  Pseudo-pure  state 
of  the  3  Carbons. 

b)  Final  state:  The  polarization  of 
the  Carbons  coupled  to  the  target 
spin  in  state  |1)  (left)  has  been 
inverted,  while  it  is  unchanged  for 
spins  coupled  to  target  spin  in  the 
|0)  state  (right). 


The  experimental  results  show  that  the  polarization  of  the  three  Carbons  is  inverted 
conditionally  on  the  state  of  the  Proton,  giving  indirect  measurement  of  its  state. 
Extending  this  scheme  to  larger  systems  with  a  much  smaller  target  to  Amplifier  spins 
ratio,  to  obtain  sensitivity  enhancement,  is  experimentally  challenging,  yet  we  can  relax 
some  of  the  requirements. 

The  signature  of  a  successful  scheme  is  that  it  produces  observable  contrast  in  a 
measurement  of  the  collective  magnetization  of  the  Amplifier.  We  define  the  contrast  as 
C  =  (Po-Pi)/  P(0),  Po  and  Pi  being  respectively  the  polarizations  for  the  target  spin  in  the 
state  |0)  or  1 1)  and  initial  polarization  P(0).  The  previous  schemes  provide  the  maximum 
contrast  C  =  2.  Allowing  for  a  lower,  but  still  observable,  contrast,  we  can  envisage 
schemes  that  do  not  require  control  over  individual  spins.  Furthermore,  it  is  not 
necessary  to  restrict  to  the  preparation  of  a  cat-state,  a  wider  class  of  entangled  states  is 
equally  useful.  Multiple  Quantum  Coherence  [3]  methods  in  NMR  offer  techniques  for 
creating  the  wanted  entangled  states. 

To  further  simplify  the  experimental  task  and  avoid  the  need  of  applying  the  inverse  of 
a  rather  complex  map,  we  can  borrow  ideas  of  Quantum  Chaos.  It  has  been  shown  that 
the  Fidelity  between  an  initial  state  evolved  under  a  unitary  chaotic  operator  and  the 
same  state  evolved  under  a  perturbed  version  of  the  same  operator  decays  exponentially 
with  the  number  of  repetitions  of  the  map,  for  strong  enough  perturbations  [4], [5].  If  we 
apply  the  perturbation  conditionally  on  the  state  of  the  target  spin,  the  two  possible 
resultant  states  will  have  an  exponentially  decreasing  overlap.  However,  the  expectation 
value  of  the  polarization  is  a  much  weaker  measurement  than  the  Fidelity  and  therefore 
we  expect  this  rate  to  simply  bound  the  contrast  growth  rate. 

In  conclusion,  we  show  that  it  is  possible  to  transfer  polarization  from  a  target  spin  to  an 
ensemble  of  spins,  through  the  creation  of  a  highly  entangled  state  with  RF  pulses  and 
the  coupling  between  the  target  spin  and  its  closest  neighbors.  The  methods  proposed 
open  the  possibility  to  a  new  class  of  devices,  where  quantum  effects,  such  as 
entanglement,  are  used  to  make  a  transition  from  microscopic  to  macroscopic 
properties. 
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When  an  ion  confined  in  an  anisotropic  bidimensional  Paul  trap  is  subjected  to  a  laser 
beam  oriented  along  an  arbitrary  direction,  the  interaction  between  its  electronic  and 
vibrational  degrees  of  freedom  is  described  by  a  time-dependent  hamiltonian  model  as  a 
consequence  of  the  lack  of  geometrical  symmetry.  Appropriately  choosing  the  laser 
frequency,  the  hamiltonian  model  turns  out  to  be  sinusoidally  oscillating  at  the  difference 
between  the  proper  frequencies  of  the  centre  of  mass  of  the  ion.  Thus,  if  the  anisotropy 
of  the  trap  is  sufficiently  small,  the  evolution  of  the  system  can  be  considered  as 
adiabatic.  In  the  context  of  this  physical  situation,  we  have  calculated  the  Berry  phase 
acquired  in  a  cycle  by  the  instantaneous  eigenstates  of  the  trapped  ion  hamiltonian. 
Suitably  choosing  both  the  initial  condition  and  a  physical  observable  we  succeed  to 
forecast  physical  effects  directly  traceable  back  to  the  accumulated  Berry  phase.  In 
particular  we  indeed  bring  to  light  that  the  mean  value  of  the  chosen  observable  after  a 
cycle  is  the  negative  of  that  calculated  at  the  same  instant  of  time  in  the  case  of  isotropic 
traps.  This  effect  demonstrates  that  and  how  the  Berry  phase  can  be  exploited  to 
evidence  the  existence  of  a  weak  anisotropy  in  a  Paul  trap. 


Amplifying  Quantum  Signals  with  a  Josephson  Bifurcation  Amplifier 
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Quantum  measurements  of  solid-state  systems,  such  as  the  readout  of 
superconducting  quantum  bits,  challenge  conventional  low-noise  amplification 
techniques.  Ideally,  the  amplifier  for  a  quantum  measurement  should  minimally  perturb 
the  measured  system  while  maintaining  sufficient  sensitivity  to  overcome  the  noise  of 
subsequent  elements  in  the  amplification  chain.  Additionally,  the  characteristic  drift  of 
materials  properties  in  solid-state  systems  necessitates  a  fast  acquisition  rate  to  permit 
measurements  in  rapid  succession.  To  meet  these  inherently  conflicting  requirements, 
we  propose  to  harness  the  sensitivity  of  a  dynamical  system  -  a  single  RF-driven 
Josephson  tunnel  junction  -  tuned  near  a  bifurcation  point. 

The  superconducting  tunnel  junction  is  the  only  electronic  dipolar  circuit 
element  with  a  non-dissipative  nonlinearity  which  remains  unchanged  at  arbitrary  low 
temperatures.  As  the  key  component  of  present  superconducting  amplifiers,  it  is  known 
to  exhibit  a  high  degree  of  stability.  Moreover,  all  available  degrees  of  freedom  in  the 
dynamical  system  participate  in  information  transfer  and  none  contribute  to  unnecessary 
dissipation  resulting  in  excess  noise. 

The  operation  of  our  Josephson  bifurcation  amplifier  is  represented 
schematically  in  Fig.  1.  The  central  element  is  a  Josephson  junction  whose  critical 
current  Io  is  modulated  by  the  input  signal  using  an  application-specific  coupling  scheme 
(input  port),  such  as  a  SQUID  loop  or  a  SSET.  The  junction  is  driven  with  an  sinusoidal 
signal  z’rf  sin(cot)  fed  from  a  transmission  line  through  a  directional  coupler  (drive  port). 
In  the  underdamped  regime,  when  the  drive  frequency  ©  is  detuned  form  the  natural 
oscillation  frequency  cop  and  when  the  drive  current  Ib'  <  z’rf  <  Ib  «  Io,  the  system  has 
two  possible  oscillation  states  which  differ  in  amplitude  and  phase  [1-2].  Starting  in  the 
lower  amplitude  state,  at  the  bifurcation  point  /rf  -  Ib  the  system  becomes  infinitely 
sensitive,  in  absence  of  thermal  and  quantum  fluctuations,  to  variations  in  Io.  The  energy 
stored  in  the  oscillation  can  always  be  made  larger  than  thermal  fluctuations  by 
increasing  the  scale  of  Io,  thus  preserving  sensitivity  at  finite  temperature.  The  reflected 
component  of  the  drive  signal,  measured  through  another  transmission  line  connected  to 
the  coupler  (output  port),  is  a  convenient  signature  of  the  junction  oscillation  state  which 
carries  with  it  information  about  the  input  signal.  The  reflected  signal  phase  <j>  is  plotted 
in  Fig  2.  as  a  function  of  the  drive  current  z’rf,  with  the  two  oscillation  states  of  the 
junction  being  labeled  “0”  and  “1”. 

This  arrangement  minimizes  the  back-action  of  the  amplifier  since  the  only 
fluctuations  felt  at  its  input  port  arise  from  the  load  impedance  of  the  follower  amplifier, 
which  is  physically  separated  from  the  junction  via  a  transmission  line  of  arbitrary 
length  and  can  therefore  be  thermalized  efficiently  to  base  temperature. 

We  present  an  experiment  that  demonstrates  the  principle  of  bifurcation 
amplification,  and  we  show  that  the  sensitivity  obtained  is  in  good  agreement  with  a 


well  established  theory  of  dynamical  transitions,  thus  making  it  relevant  for  stringent 
tests  of  superconducting  qubits  and  gates. 


iRFsin(cot+$rt0iRF)) 


iRF  sin(<ot) 


Figure  1:  Schematic  diagram  of  the  Josephson  bifurcation 
amplifier.  A  junction  with  critical  current  I0  coupled  to  the 
input  port  is  driven  by  an  RF  signal.  In  the  vicinity  of  the 
dynamical  bifurcation  point,  the  phase  of  the  resulting 
oscillation  depends  critically  on  the  input  signal  and 
manifests  itself  in  the  reflected  signal  phase  <j>. 


drive  current  iRF  / 10 


Figure  2:  Reflected  signal  phase  <f>  as  a  function  of  drive  current 
/rf/Io-  Io=1.172pA;  (<»p-o))/co=0.122.  Symbols  denote  the  mode 
of  <j>  with  up  and  down  triangles  corresponding  to  increasing 
and  decreasing  current,  respectively.  The  dotted  line  is 
The  calculated  bifurcation  points,  IB-  and  IB,  are  marked  on  the 
horizontal  axis.  The  0  and  1  phase  states  are  reminiscent  of  the 
superconducting  and  dissipative  states  of  the  DC  current  biased 
junction. 
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Superconducting  circuits  containing  Josephson  junctions  are  very  promising  candidates 
for  the  implementation  of  solid-state  quantum  bits  or  qubits  [1].  Complex  single-qubit 
operations  have  already  been  reported  [2,3],  as  well  as  coherent  dynamics  of  two  coupled 
qubits  [4]  and  realization  of  a  quantum  gate  [5].  Coupling  a  qubit  to  a  harmonic  oscillator 
is  interesting  from  a  fundamental  point  of  view  to  generate  and  study  non-classical  states 
of  the  oscillator,  and  is  relevant  in  quantum  information  as  well,  since  coupling  many 
qubits  via  a  harmonic  oscillator  has  been  proposed  [6],  Here  we  present  recent 
measurements  demonstrating  the  coupling  of  a  flux-qubit  to  the  plasma  mode  of  the  DC 
Squid  which  at  the  same  time  is  used  to  measure  the  qubit  state.  This  coupling  is 
manifested  by  the  appearance  of  two  side-band  resonances  around  the  bare  qubit  peak. 
By  performing  two-pulse  experiments,  we  demonstrate  that  these  additional  resonances 
are  indeed  excitations  of  the  coupled  system.  We  also  observe  Rabi  oscillations  between 
the  coupled  |qubit,Squid>  states,  thus  demonstrating  entanglement  between  the  states  of 
two  superconducting  circuits.  We  use  the  qubit  to  measure  intrinsic  properties  of  the 
plasma  mode  :  temperature,  relaxation  time.  Conversely,  we  also  measure  the  qubit's  state 
via  the  plasma  mode.  These  results  indicate  that  complex  manipulation  of  entangled 
states  similar  to  cavity  QED  or  trapped-ions  experiments  is  within  reach  with 
superconducting  circuits. 
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A  superconducting  wire  can  be  considered  as  quasi-one  dimensional  (ID)  if  its 
characteristic  transverse  dimension  Vo  (o  being  the  cross  section)  is  smaller  than  the 
coherence  length  ^(T).  The  shape  of  the  bottom  part  of  the  resistive  transition  R(T)  of  a  ID 
superconducting  strip  is  described  by  the  model  of  phase  slips  activation.  If  the  wire  is 
infinitely  long,  then  there  is  always  a  finite  probability  that  a  ‘small’  part  of  the  sample  is 
instantly  driven  normal:  locally  the  modulus  of  the  order  parameter  |Aj  goes  to  zero  and  its 
phase  cp  ‘slips’  by  2 n. 

Free  energy  of  a  current-carrying  ID  superconductor  is  described  by  so  called  ’tilted 
washboard’  potential  (Fig.  1).  The  black  circle  denotes  a  thermodynamically  metastable 
state  of  a  macroscopically  coherent  superconductor.  The  system  can  relax  to  a  state  of 
lower  energy.  There  are  two  possible  ways  for  a  superconductor  to  decrease  its  energy. 
First,  due  to  thermal  excitation  it  can  ‘jump’  over  the  energy  barrier  AF  to  a  new  potential 
minimum  [1]:  solid  arrow  in  Fig.  1.  Second,  the  system  can  tunnel  through  the  energy 
barrier  [2]:  dotted  arrow  in  Fig.  1.  In  the  limit  of  infinitely  small  supercurrent  the  height  of 
the  barrier  AF  is  the  energy  difference  between  normal  and  superconducting  states  of  a 
system  of  the  smallest  possible  size  Q  =  £  ■  o:  AF  =  Bc2  •  Q,  Bc  being  the  critical  magnetic 
field.  Corresponding  decrease  of  energy  is  accompanied  by  the  change  in  phase  by  2n  (Fig. 
1).  It  can  be  shown  [1],  that  in  case  of  thermally  activated  phase  slips  the  effective 
resistance  R(T)  ~  exp  (-AF  /  kBT).  For  a  not  too  long  wire  L  >  i|(T)  (only  single  phase  slip 
event  can  happen  at  a  time)  with  normal  state  resistance  Rn  the  quantum  tunneling 
contribution  is  [2]:  R(T)  /  RN  =  (£,  /  L)  •  (to  •  Tqps),  where  To  ~  h  /  A  is  the  duration  of  a 
phase  slippage,  and  Fqps  -  B  exp  (-Sqps)  is  the  rate  of  quantum  phase  slip  events,  where  B 
~  (Sqps  /  To)  •  (L  / 1),  SQPS  =  A-(Rq  /  RN)-(L  /  A  ~  1,  Rq  =  h  /  4e2  =  6.47  kQ. 

Fig.  1.  Dependence  of  the  free 
energy  F  vs.  superconducting  phase 
tp  of  a  ID  current-carrying 
superconductor.  The  system 
(represented  by  a  circle)  can 
change  its  quantum  state  in  two 
ways:  by  thermally  activated  phase 
slips  (solid  arrow),  or  by  quantum 
tunneling  (dotted  arrow).  T  stands 
for  the  rates  of  the  corresponding 
processes. 


Several  experimental  attempts  have  been  made  to  detect  the  QPS  mechanism  in  narrow 
superconducting  channels:  ultra-thin  In  and  Pb-In  wires  [3]  and  Mo-Ge  film  evaporated  on 
top  of  a  carbon  nanotube  [4].  The  results  are  very  intriguing,  but  the  matter  is  far  from 
being  settled.  The  motivation  of  this  work  was  to  study  the  manifestation  of  the  QPS 
mechanism  in  the  same  superconducting  wire  as  a  function  of  its  cross  section. 

We  have  developed  a  method  of  progressive  reduction  of  transverse  dimensions  of  e- 
beam  lift-off  fabricated  nanostructures  by  ion-beam  sputtering  [5].  The  method  enables 
galvanomagnetic  measurements  of  the  same  sample  in  between  the  sessions  of  etching. 
Sputtering  can  be  considered  as  an  erosion  of  the  surface  due  to  bombardment  of  primary 
ions.  The  method  is  very  promising,  for  it  lets  us  directly  follow  changes  in 
superconductive  transition  in  a  ID  superconductor  along  with  successive  reduction  of  its 
cross  section.  AFM  control  revealed  that  sputtering  makes  the  sample  surface  smoother  and 
more  homogenous,  eliminating  inevitable  imperfections  while  original  lift-off  process.  For 
not  too  narrow  Al  wires  the  shape  of  R(T)  transition  can  be  perfectly  fitted  by  the  thermal 
activation  model  [1].  While  for  the  wires  that  after  few  sputtering  sessions  have  reached  the 
effective  diameter  Va  <  35  nm  a  low  temperature  ‘foot’  develops,  which  can  be  associated 
with  the  QPS  phenomena  (Fig.  2). 

This  work  has  been  supported  by  the  EU  Commission  "ULTRA- ID"  project  NMP4- 
CT-2003-505457  "Experimental  and  theoretical  investigation  of  electron  transport  in  ultra¬ 
narrow  1 -dimensional  nanostructures". 
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Fig.  2.  Resistive  transition  R(T)  of  a  5  pm  long  same  Al  wire  before  sputtering  (left  curve) 
and  after  (right).  High-temperature  part  of  the  transition  can  be  described  by  the  thermal  PS 
activation  mechanism  (solid  lines),  while  for  the  thinner  wire  the  bottom  ‘foot’  can  be  fitted 
by  the  QPS  model. 
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In  this  paper  we  study  a  Josephson  junction  (JJ)  (kept  at  cryo-temperature  T)  when 
connected  to  an  external  current  source  I(t),  via  a  twisted  pair  with  impedance  Ztw.  See 
Fig  1(a).  The  generic  features  of  such  systems  have  been  studied  extensively  by  e.g. 
Kautz  and  Martinis  [1],  Here  however  we  are  interested  in  the  quantitative  behaviour  of 
our  specific  system.  We  study  the  probability  distribution  P  of  switching  of  the  JJ,  as  a 
function  of  the  current  pulse  from  the  source.  Due  to  quantum  effects  and  various  noise 
sources,  such  as  resistive  leads  at  finite  temperatures,  this  distribution  will  have  a  finite 
width.  In  a  qubit  system,  such  as  the  Quantronium  [2],  it  is  necessary  that  this  width  is 
small  in  order  to  clearly  distinguish  between  different  quantum  states,  and  eventually 
realise  the  single  shot  measurement.  Here  we  will  compare  simulated  and  measured 
values  of  this  width. 

To  numerically  calculate  the  function  P,  we:  (i)  model  the  junction  with  an  ideal 
element  in  parallell  with  a  capacitor  Cj.  (ii)  for  the  current  pulse  we  use  a  technique 
referred  to  as  pulse  and  hold  [3],  with  a  fast  switching  pulse  (with  magnitude  Ip  and 
duration  xp)  followed  by  a  longer  hold  pulse,  (iii)  We  model  Ztw  with  a  resistor  R]  in 
parallell  with  another  resistor  R2  and  a  capacitor  C2  (the  impedance  corresponding  to 
this  setup  is  called  ZRRC,  see  Fig  1(b).  Zrrc  is  fitted  to  match  Ztw  as  closely  as  possible 
(see  Fig.  2(a),  (b)  for  the  real  and  imaginary  part  of  the  measured  impedance  (solid  line) 
and  the  model  ZRRC  (dashed  line),  (iv)  The  thermal  noise  from  Ztw  we  model  with  the 
noise  currents  Ini>2  =  NvTn/R12  -A  where  N  is  a  Gaussian  distributed  random  number 
and  A  is  a  constant. 

Applying  Kirchoffs  rules  and  the  Josephson  relations  to  the  circuit  in  Fig  1(a)  with  the 
model  impedance  in  Fig  1(b)  we  find  a  third  order  differential  equation  for  the  phase 
over  the  junction.  This  description  is  valid  here  since  Ej »  Ec  (Ec~0.5  K,  Ej~6.3  K). 

The  width  AI  of  P  we  define  as  AT=a"’(r),  where  r=a(AI)=|P(Isw+AI/2)-P(Isw-AI/2)|  is  the 
discriminating  powe  and  Isw  is  the  switching  current,  P(Isw)=l/2. 

In  Fig.  3(a)  we  plot  AI/ISW,  and  in  Fig.  3(b)  we  plot  AI  as  a  function  of  Xp.  The  stars 
connected  with  solid  lines  are  measured  data  with  0.025=T=0.7  K  (we  do  not  bother  of 
identifying  the  temperature  for  the  curves),  and  the  stars  connected  with  dashed  lines 
are  simulated  data  with  Tn=0.03,  0.7  and  4.0  K. 

For  Xp=10  ps,  the  simulated  data  agrees  well  with  the  measured  ones  if  Tn  =  0.03  K.  For 
Xp=0.1  ps  the  AI/ISW  values  agree  for  Tn=4  K  whereas  the  AI  values  do  not  agree  at  all. 
Increasing  Tn  further  will  (most  likely)  not  increase  AI  in  any  considerable  way.  A  more 


thorough  analysis  indicates  that  the  big  difference  in  the  AI  values  for  Xp=0. 1  ps  is  at 
least  not  due  to  the  simplification  of  the  impedance. 

(a)  (b) 
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Several  quantum-control  experiments  based  on  Josephson  junctions  measure  the 
switching  current  to  determine  the  quantum  state  of  the  circuit  [1][2],  Some  of  these 
experiments  rely  on  a  pulsed  readout  where  a  rectangular  current  pulse  is  applied  to  the 
circuit  and  the  voltage  is  measured  to  determine  switching.  The  sample  is  typically 
mounted  in  a  cryostat  and  small  currents  have  to  charge  the  capacitance  of  the 
measurement  leads,  delaying  the  signal  at  the  readout  stage.  Hence,  a  switching  in  the 
very  last  instance  of  the  pulse  cannot  build  up  enough  voltage  needed  to  trigger  a 
detector.  Here  we  present  a  method  to  circumvent  this  problem.  A  pulse  of  the  order  of 
the  switching  current,  called  the  switch  pulse,  is  immediately  followed  by  a  hold  level 
of  smaller  amplitude.  Its  purpose  is  to  latch  the  state  of  the  circuit  for  the  time  that  is 
needed  until  the  readout  is  finished.  With  this  hold  level,  a  late  switch  can  be  detected 
and  switching  ambiguities  are  excluded. 

After  switching,  the  new  stable  operating  point  of  the  circuit  is  defined  by  the 
intersection  of  the  sample  IV  and  the  loadline,  set  by  the  bias  resistor  The  choice  of 

the  bias  resistor  depends  on  two  things.  On  the  one  hand,  the  operating  point  should  be 
at  a  voltage  as  low  as  possible,  yet  large  enough  to  distinguish  switch  from  no-switch. 
In  this  case,  the  power  dissipated  in  the  junction  is  minimal  since  a  current  through  the 
device  is  minimal.  On  the  other  hand,  a  latching  behaviour  is  only  obtained  if  the  slope 
of  the  loadline,  which  is  -l/R^,  is  smaller  than  the  negative  differential  resistance  of  the 

junction  IV. 

The  experimental  setup  is  shown  in  figure  1.  An  arbitrary  waveform  generator 
(AWG)  together  with  an  attenuator  and  a  bias  resistor  form  the  bias  circuit.  It  is 
connected  to  the  sample  via  twisted  pair  cables.  High  frequency  noise  is  filtered  at  room 
temperature  with  a  15MHz  low-pass  MiniCircuit  filter.  The  sample  is  placed  in  a 
dilution  refrigerator  with  a  base  temperature  of  25mK.  The  detection  of  the  response 
voltage  from  the  sample  is  done  after  the  bias  resistor.  While  the  sample  stays  in  the 
superconducting  state,  the  main  portion  of  the  applied  voltage  is  dropped  across  and 

only  a  small  voltage  is  registered.  If  the  sample  switches  to  the  normal  state  with  a  high 
quasi-particle  resistance  (much  larger  than  the  bias  resistor),  the  main  voltage  drop  is 
across  the  sample  and  a  larger  voltage  than  in  the  no-switch  case  is  detected.  An  event 
counter  with  a  threshold  voltage  set  between  the  switch  and  no-switch  signal  counts  the 
switching  events.  Since  the  minimum  threshold  that  can  be  programmed  into  the 
counter  is  50m V,  the  signals,  which  are  below  the  gap  of  about  400pV  have  to  be 
amplified.  A  home-built  pre-amplifier  with  a  gain  of  500  is  responsible  for  this  task. 
Finally,  a  fast  digital  sampling  oscilloscope  (DSO)  with  large  memory  can  record 
complete  measurement  traces  of  10^  pulses  and  therefore  offers  the  possibility  of  doing 
correlation  analysis. 

The  experiments  were  carried  out  on  SQUID  junctions,  see  inset  of  fig.  1.  The 
SQUID  allows  us  to  perform  measurements  in  a  wide  range  of  switching  currents,  from 
the  maximum  switching  current  to  zero,  by  applying  a  magnetic  field  to  the  SQUID. 


Thus,  the  ratio  of  the  two  important  energies  of  this  system,  which  are  the  Josephson 
energy  Ej  and  the  charging  energy  E^,  can  be  continuously  varied. 

In  a  typical  experiment  a  total  of  Nf=104  identical  pulses,  each  separated  by  a 
wait  time  xw,  are  applied  to  the  sample.  The  hold  level  is  set  that  it  never  switches  the 

sample.  This  level  is  maintained  throughout  the  whole  experiment.  The  amplitude  of  the 
switch  pulse  is  initially  at  a  level  that  no  pulse  causes  a  switch.  It  is  then  continuously 
increased  until  all  Nt  pulses  switch  the  sample.  The  switching  probability  is  obtained 

from  P=N/Nt,  where  N  is  the  number  of  switches.  From  the  width  of  this  switching 

distribution,  the  resolution  AI  of  such  a  detector  is  obtained.  We  study  the  normalized 
resolution,  AI/I,  where  I  is  the  current  value  at  P=0.5.  For  small  switching  currents  of 
the  order  of  12nA,  we  find  a  resolution  of  less  than  10%  for  switching  pulses  of  100ns 
duration.  For  longer  pulses  or  pulses  with  higher  switching  currents,  the  resolution 
decreases  even  further. 

After  each  switch,  the  system  needs  time  to  relax  to  its  initial  state.  Therefore, 
the  wait  time  has  to  be  chosen  long  enough  in  order  that  a  pulse  won’t  influence  the 
outcome  of  a  following  pulse.  We  made  a  correlation  analysis  where  we  calculate  the 
auto-correlation  of  a  pulse-sequence  as  a  function  of  the  wait  time,  while  the  switching 
probability  is  held  constant  at  P=0.5,  see  fig.  2.  We  find  a  substantial  auto-correlation 
for  wait  times  shorter  than  40ps.  A  long  decay  time  of  the  auto-correlation  function 
might  be  explained  by  excessive  heating  due  to  quasi  particle  tunnelling  in  the  junction 
when  it  switches  above  the  gap  voltage.  However,  since  in  our  case  we  keep  the  voltage 
below  the  gap,  only  small  currents  discharge  the  capacitance  of  the  leads.  Therefore  it  is 
possible  that  the  junction  does  not  always  retrap  before  a  new  pulse  is  applied. 
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We  are  conducting  an  experimental  study  on  the  effect  of  high  impedance 
environment  on  a  Cooper  pair  transistor  (CPT).  The  CPT  consists  of  two  small- 
capacitance  Josephson  tunnel  junctions  in  series  with  a  gate  electrode  coupled  through 
a  capacitance  Cg  to  a  central  island.  In  small-capacitance  CPT  the  charging  energy 
Ece  =  e2/  2Cs,  where  Cs  ~  2Cj  +  Cg  is  the  total  capacitance  of  the  island,  becomes 
relevant  at  low  temperature,  and  charging  effects  influence  the  transport  properties.  In 
order  to  observe  Coulomb  blockade  in  a  single  Josephson  junction  and  also  in  a  CPT, 
care  must  be  taken  to  isolate  the  junctions  from  their  high  frequency  electromagnetic 
environment.  The  classical  Josephson  effect,  realized  in  the  low  impedance 
environment  (Re[Z(©)]  <  Rq  =  h/4e2  -6.45  kQ),  is  suppressed  in  the  high  impedance 
environment  where  there  are  large  quantum  fluctuations  of  the  phase.  In  a  high 
impedance  environment,  the  current  voltage  (I-V)  curve  of  the  CPT  does  not  present  a 
supercurrent  feature,  but  a  Coulomb  blockade  feature,  and  a  region  of  negative 
differential  resistance  which  is  the  signature  of  coherent  tunnelling  of  Cooper  pairs 
[1],  In  order  to  observe  the  Coulomb  blockade  of  Cooper  pair  in  a  CPT,  the  main 
requirement  on  the  environment  is  that  its  impedance  (Re[Z(ra)])  shunting  the  CPT 
must  be  much  larger  than  the  quantum  resistance  Rq.  In  previous  works  [2],  the 
impedance  of  the  leads  was  increased  by  using  small  chromium  resistors  with  a 
resistance  of  50  kQ,  located  on  the  chip  close  to  the  junctions.  In  our  experiments,  as 
in  those  of  Watanabe  [3],  the  high  impedance  of  the  leads  is  obtained  by  biasing  the 
CPT  with  four  SQUID  arrays.  The  SQUID  geometry  allows  the  tuning  in  situ  of  the 
leads  impedance  over  several  orders  of  magnitude  (104  Q<Re[Z(©=0)]<108  Q), 
providing  different  electromagnetic  environments  to  the  CPT.  Note  that  the  junctions 
defining  the  CPT  do  not  have  a  SQUID  geometry,  and  thus  are  unaffected  by  the 
relatively  small  magnetic  fields  (the  flux  O=O0/2  at  5  mT,  where  O0=h/2e  =2x1  O'1 5 
Wb  is  the  flux  quantum)  needed  to  tune  the  SQUIDs. 

Measurements  were  performed  in  a  dilution  refrigerator  with  a  base  temperature  of 
15mK.  No  low  temperature  filtering  is  implemented  in  the  cryostat.  However,  when 
the  arrays  are  in  the  high  impedance  state,  they  act  as  good  filters  for  electromagnetic 
fluctuations  guided  to  the  sample  via  the  leads.  The  I-V  characteristics  of  the  CPT  are 
measured  in  a  four  point  configuration  (Fig.  1).  One  pair  of  SQUID  leads  are  used  to 
apply  a  symmetric  bias  and  measure  the  current,  the  other  pair  to  measure  the  voltage 
across  the  CPT.  Each  array  has  70  SQUIDs,  with  an  area  of  0.03  pm2  per  junction. 
The  CPT  junctions  are  designed  to  have  an  area  of  0.01  pm2  each.  CPTs  with  different 
ratio  Ej/Eci  have  been  measured. 

In  several  experiments  performed  in  a  low  impedance  environment  [4,5],  it  has 
been  shown  how  the  2e-periodicity  of  the  supercurrent  in  the  gate  induced  charge  is 
suppressed  by  non-equilibrium  quasiparticles.  In  order  to  study  this  effect  in  a  high 
impedance  environment,  we  fixed  the  bias  voltage  across  the  SQUID  arrays  at  a  value 
of  50  pV  and  measured  the  current  and  voltage  of  the  transistor  as  function  of  Vg.  The 


tunneling  of  charges  in  and  out  the  island  results  in  an  oscillating  behaviour  of  both 
the  current  and  voltage  across  the  transistor  (Fig  2).  It  is  interesting  to  notice  how  the 
character  of  the  oscillations  changes  as  the  impedance  of  the  environment  is 
increased.  In  fact,  as  the  magnetic  field  reduces  the  coupling  energy  of  the  SQUID 
arrays  (increasing  their  impedance),  the  Vg  dependence  changes  from  e-periodic  to 
2e-periodic.  The  change  in  periodicity  is  observed  in  both  current  and  voltage 
oscillations.  This  suggests  that  non-equilibrium  quasiparticles  present  in  the  system 
and  responsible  for  the  e-periodicity,  can  be  suppressed  by  high  impedance  leads. 


Figure  1 


Figure  2 
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Fig  1:  A  SEM  micrograph  of  the  CPT  biased  by  4  SQUID  arrays  and  measurement 
scheme. 

Fig  2:  Voltage  and  current  modulation  as  function  of  Vg  for  different  Re[Z(co)].  From 
top  to  bottom:  55  kQ,  450  kO,  10  MQ,  20  Mfi.  For  this  CPT  Ej/Ecr  =0.62. 
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Quantum  dots  for  single  photon  and  photon  pair  technology. 
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Quantum  dots  have  emerged  as  a  popular  platform  for  developing  single  photon 
emitting  devices,  an  essential  component  in  many  areas  of  quantum  information 
technology.  Much  progress  has  already  been  made,  with  the  separate  realisation  of 
electrically  driven  [1],  and  quite  efficient  [2]  devices.  These  demonstrations  rely  on  the 
unique  Coulomb  energy  for  each  exciton  complex.  Usually  the  radiative  decay  of  the 
neutral  single  exciton  is  spectrally  filtered  to  output  only  a  single  photon  per  cycle. 

Further  consideration  raises  an  interesting  question;  what  of  the  other  exciton 
complexes  formed  during  the  decay  of  the  initially  excited  exciton  state?  The  initial  state 
is  not  unique,  due  to  the  probabilistic  nature  of  carrier  capture,  and  thus  radiative  decay 
can  occur  via  a  multi  photon  cascade.  It  is  also  possible  that  non-radiative  (dark)  initial 
exciton  states  can  be  formed  when  the  spins  of  the  component  electron  and  heavy  hole 
have  the  same  sign.  These  processes  clearly  must  have  an  impact  on  operation,  and 
device  design  must  embrace  any  advantages  that  exist.  Thus,  we  present  time  integrated 
and  time  resolved  micro-photoluminescence  experiments  that  reveal  the  broader 
properties  of  single  photon  and  multi-photon  emission  from  InAs  quantum  dots. 

The  dots  studied  here  are  fabricated  by  growing  a  single  layer  of  InAs  at  the 
critical  thickness  for  island  formation  in  a  GaAs  matrix,  so  that  an  areal  density  of  less 
than  lp"2  is  achieved.  These  dots  are  thought  to  be  rather  small  in  size,  and  show 
relatively  few  emission  lines  in  the  spectra,  due  to  the  lack  of  a  confined  excited  exciton 
state.  We  characterise  the  radiative  decay  of  the  biexciton  state,  and  show  reduced  jitter 
for  single  photon  emission  of  biexciton  photons,  and  correlation  with  the  linear 
polarisation  of  the  subsequent  exciton  photon.  The  correlation  is  shown  to  exist  due  to 
the  lifting  of  the  degeneracy  of  the  intermediate  optically  active  exciton  level,  attributed 
to  in-plane  asymmetry  of  the  electron  and  hole  wavefunctions.  The  choice  of 
polarisation  for  each  linearly  polarised  pair  of  photons,  is  shown  to  be  a  random  for  each 
emission  cycle,  and  we  discuss  how  this  might  be  used  to  implement  optical  quantum 
key  generation  [3]. 

We  also  observe  evidence  for  non-radiative  processes,  including  scattering 
between  optically  active  exciton  sub-levels,  and  the  formation  of  very  long  lived 
optically  inactive  states.  The  latter  is  inferred  from  the  lifetimes  of  predominantly 
optically  inactive  mixed  states  [4],  created  by  in-plane  magnetic  fields.  This  is  important 
for  single  photon  emitters,  where  dark  states  can  block  emission  for  whole  periods,  and 
increase  jitter  by  subsequent  excitation  into  higher  order  exciton  states. 
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Phonon-induced  decay  of  the  electron  spin  in  quantum  dots 
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Phase  coherence  of  spin  in  quantum  dots  is  of  central  importance  for  spin-based 
quantum  computation  in  the  solid  state  [1,2].  Sufficiently  long  coherence  times  are 
needed  for  implementing  quantum  algorithms  and  error  correction  schemes.  If  the  qubit 
is  operated  as  a  classical  bit,  its  decay  time  is  given  by  the  spin  relaxation  time  T\, 
which  is  the  time  of  a  spin-flip  process.  For  quantum  computation,  however,  the  spin 
decoherence  time  T2  —  the  lifetime  of  a  coherent  superposition  of  spin  up  and  spin 
down  states  —  must  be  sufficiently  long.  In  semiconductor  quantum  dots,  the  spin 
coherence  is  limited  by  the  dot  intrinsic  degrees  of  freedom,  such  as  phonons,  spins  of 
nuclei,  excitations  on  the  Fermi  surface  ( e.g .  in  metallic  gates),  fluctuating  impurity 
states  nearby  the  dot,  electromagnetic  fields,  etc.  It  is  well  known  (and  experimentally 
verified)  that  the  T\  time  of  spin  in  quantum  dots  is  extremely  long,  extending  up  to  100 
ps.  The  decoherence  time  7),  in  its  turn,  is  limited  by  both  spin-flip  and  dephasing 
processes,  and  can  be  much  smaller  than  7),  although  its  upper  bound  is  T2<  27). 
Knowledge  of  the  mechanisms  of  spin  relaxation  and  decoherence  in  quantum  dots  can 
allow  one  to  find  regimes  with  the  least  spin  decay. 

In  our  work  [3],  we  study  spin  relaxation  and  decoherence  in  a  GaAs  quantum  dot  due 
to  spin-orbit  interaction.  We  derive  an  effective  Hamiltonian  which  couples  the  electron 
spin  to  phonons  or  any  other  fluctuation  of  the  dot  potential.  We  show  that  the  spin 
decoherence  time  T2  is  as  large  as  the  spin  relaxation  time  7),  under  realistic  conditions. 
For  the  Dresselhaus  and  Rashba  spin-orbit  couplings,  we  find  that,  in  leading  order,  the 
effective  magnetic  field  can  have  only  fluctuations  transverse  to  the  applied  magnetic 
field.  As  a  result,  T2=2T\  for  arbitrarily  large  Zeeman  splittings,  in  contrast  to  the 
naively  expected  case  T2«T\.  We  show  that  the  spin  decay  is  drastically  suppressed  for 
certain  magnetic  field  directions  and  values  of  the  Rashba  coupling  constant.  Finally, 
for  the  spin  coupling  to  acoustic  phonons,  we  show  that  T2=2T\  for  all  spin-orbit 
mechanisms  in  leading  order  in  the  electron-phonon  interaction. 
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The  spin  of  a  single  electron,  confined  in  a  semiconductor  quantum  dot,  provides  a  natural 
two-level  system  suitable  as  a  qubit  in  a  quantum  computer  [1],  We  demonstrate  here 
recent  experimental  progress  towards  realization  of  an  electron-spin  based  quantum 
computer:  isolation  of  a  single  electron,  identification  of  the  two-level  spin  system  and 
single-shot  read-out  of  the  spin  state  of  a  single  electron. 

Our  device,  defined  electrostatically  in  a  GaAs/AlGaAs  heterostructure,  consists  of  two 
coupled  quantum  dots  in  close  proximity  to  a  quantum  point  contact  (QPC).  Using  this 
QPC  as  a  charge  detector,  we  show  that  we  can  deplete  the  two  dots  to  contain  just  a  single 
electron  each  [2],  Even  in  this  few-electron  regime,  the  resonant  current  through  the  two 
dots  in  series  is  still  measurable. 

By  applying  an  in-plane  magnetic  field  and  performing  large-bias  tunneling  spectroscopy, 
we  can  directly  measure  the  Zeeman  splitting  of  the  orbital  ground  state,  thereby 
identifying  the  spin  qubit  two-level  system  [3],  Furthermore,  we  show  that  the  quantum  dot 
can  be  operated  as  an  electrically  tunable  spin  filter.  The  polarization  of  the  filter  is 
opposite  for  electrons  tunneling  to  the  one-  and  two-electron  ground  state,  and  is  estimated 
to  be  very  close  to  100%  in  magnitude  [4], 

To  use  this  system  as  a  spin  qubit,  an  essential  ingredient  is  a  read-out  device,  i.e.  a  way  to 
measure  the  orientation  of  a  single  electron  spin.  We  use  the  Zeeman  energy  difference 
between  the  two  spin  orientations  to  achieve  spin-to-charge  conversion;  spin  “up”  has  a 
low  energy  and  is  trapped  on  the  dot,  whereas  spin  “down”  has  a  higher  energy  and  can 
escape  to  the  leads.  By  measuring  the  charge  dynamics,  using  the  nearby  QPC  with  a 
bandwidth  of  -100  kHz,  we  are  able  to  perform  single-shot  read-out  of  the  spin.  We  use 
this  method  to  determine  the  relaxation  time  of  the  spin  qubit,  and  find  T\  =  0.85  ms  for  an 
in-plane  magnetic  field  of  8  T. 
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Research  of  low  dimensional  organic  conductors,  polymers  and  magnetic  metallic 
nanostructures  has  a  long  tradition  in  our  institute.  We  also  perform  theoretical  studies  of  such 
systems.  Recently  our  interests  are  focused  on  electronic  transport  through  nanostructures  and 
single  molecules.  In  these  systems  electronic  correlations  and  interference  processes  play  an 
important  role,  leading  in  some  cases  to  the  Kondo  and  the  Fano  resonance.  Magnetic  systems 
are  every  promising  for  their  potential  applications  in  electronics.  Interplay  between  the 
charge  and  the  spin  degree  of  freedom  is  a  relevant  aspect  in  our  research.  These  problems 
will  be  presented  for  the  coherent  transport  through  simple  models  of  molecules  (for  a  single 
atom,  a  two-atomic  molecule  and  a  carbon  nanotube).  Magnetoresistance  is  an  important 
quantity  for  operating  magnetoelectronics  devices,  therefore  its  features  will  be  presented  in 
detail.  The  magnetoresistance  is  positive  for  uncorrelated  electron  transport;  however, 
correlations  can  influence  the  magnetoresistance  changing  its  value  and  even  its  sign.  The 
problem  of  correlations  and  many-body  effects  will  be  also  discussed  for  strong  non¬ 
equilibrium  electronic  transport  through  molecules. 
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In  this  paper  we  consider  the  behaviour  of  a  quantum  regime  SQUID  ring  (here,  one 
Josephson  weak  link  enclosed  by  a  thick  superconducting  ring)  from  the  viewpoint  of 
generating  superpositions  of  macroscopically  distinct  states  and  in  the  squeezing  of 
coherent  states  [1,2].  We  show  that  because  of  the  Josephson  coupling  term  in  the 
SQUID  ring  Hamiltonian  both  are  possible.  On  this  basis  we  consider  that  the  SQUID 
ring  constitutes  a  strong  contender  for  application  in  future  quantum  circuit 
technologies.  In  order  to  make  our  calculations  correspond  to  a  more  realistic  situation 
we  also  consider  the  effect  on  the  SQUID  ring  of  a  dissipative  (decohering)  external 
environment,  taking  the  standard  open  systems  approach  commonly  used  in  the  field  of 
Quantum  Optics.  We  show  that  even  with  such  a  decohering  environment  in  place  it  is 
still  possible  to  utilise  SQUID  ring  for  creating  superpositions  and  squeezing  provided 
the  dissipation  is  not  too  great. 

As  another  example  of  applying  the  techniques  familiar  to  the  field  of  quantum  optics, 
we  consider  the  problem  of  entangling  a  SQUID  ring  coupled  inductively  to  an 
electromagnetic  field  oscillator  mode.  We  demonstrate  by  solving  the  time  dependent 
Schrddinger  equation  that  we  can  create  a  state  of  maximal  entanglement  between  them. 
We  show  that  this  state  can  be  generated  and  controlled  through  the  application  of  a 
time  dependent  external  magnetic  flux  bias  applied  to  the  ring.  When  this  flux  is 
adiabatically  ramped  through  an  avoided  crossing  we  are  able  to  create  a  reasonably 
good  approximation  to  a  maximally  entangled  state.  However,  we  find  that  this  state  of 
entanglement  can  be  enhanced  considerably  by  using  a  non-adiabatic  flux  pulse.  We 
extend  these  calculations  to  include  the  effects  of  dissipative  environments  coupled  to 
the  SQUID  ring.  We  demonstrate  that  even  including  such  decohering  environments  it 
is  possible  to  maintain  a  sufficient  degree  of  entanglement  over  time  scales  that  would 
be  required  for  quantum  technologies. 
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In  this  work,  and  by  analogy  with  Quantum  Optics,  we  discuss  the  highly  non- 
perturbative  properties  of  SQUID  rings  (here,  a  Josephson  weak  link  enclosed  by  a 
thick  superconducting  ring)  in  the  quantum  regime.  In  particular,  we  investigate  the 
possibility  of  frequency  down  conversion,  via  a  quantum  mechanical  SQUID  ring, 
between  quantised  electromagnetic  field  oscillator  modes.  We  show  that  because  of 
these  non-perturbative  properties  (essentially  due  to  the  cosinusoidal  Josephson 
coupling  term  in  the  SQUID  ring  potential)  that  multi-photon  down  conversion  can  take 
place  between  an  input  photon  and  output  photons  at  sub-integers  of  the  input 
frequency.  In  addition,  these  output  photons  are  created  in  an  entangled  state  with  the 
degree  of  entanglement  adjustable  through  the  static  magnetic  flux  applied  to  the 
SQUID  ring.  Given  that  it  appears  from  these  calculations  that  the  number  of  output 
photon  states  can  be  large,  we  consider  the  application  of  such  entangled  photons  in 
quantum  technologies. 
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SQUID  rings  irradiated  with  nonclassical  electromagnetic  fields,  are  studied.  The 
Josephson  current  is,  in  this  case,  a  quantum  mechanical  operator,  whose  expectation 
value  with  respect  to  the  density  matrix  describing  the  external  electromagnetic  field 
yields  the  experimentally  observed  current. 

Two  mesoscopic  SQUID  rings  A  and  B,  which  are  far  from  each  other,  are  considered. 
A  source  of  two-mode  nonclassical  microwaves  irradiates  the  two  rings  with  correlated 
photons.  Classically  correlated  (separable)  and  quantum  mechanically  correlated 
(entangled)  microwaves  are  considered,  and  their  effect  on  the  Josephson  currents  is 
quantified,  using  the  external  field  approximation.  Results  for  different  examples  that 
involve  microwaves  in  number  states  and  coherent  states  are  derived.  It  is  shown  that 
the  quantum  currents  and  the  quantum  statistics  of  the  electron  pairs  tunneling  through 
the  distant  Josephson  junctions  in  the  two  rings  are  correlated. 
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A  system  comprising  of  two  squid  rings  irradiated  with  entangled  electromagnetic 
fields  is  considered.  The  Hamiltonian  describing  the  system  is: 

H=Jaz  +y0'z  +^)  +  a2(a;a2  +  ^)  +  AAaf(al  +  af)  +  ABaf  (a2  +af) 

The  time  evolution  of  this  system  is  studied  numerically.  Both  separable  and 
entangled  electromagnetic  fields  are  considered  and  their  effects  on  the  currents  is 
calculated.  It  is  shown  that  the  currents  in  the  two  Squid  rings  are  correlated. 
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We  implement  the  impedance  measurement  technique  (IMT)  for  characterization  of 
interferometer-type  superconducting  qubits.  In  the  framework  of  this  method,  the  interferometer 
loop  is  inductively  coupled  to  a  superconducting  high-quality  tank  circuit.  Conclusive 
information  about  qubits  is  obtained  from  the  read-out  of  the  tank  properties.  We  show  that  the 
IMT  is  a  powerful  tool  to  study  quantum  properties  of  superconducting  structures. 
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The  phase  diagram  of  a  2D  Josephson  junction  array  with  large  substrate  resistance, 
described  by  a  quantum  XY  model,  is  studied  by  means  of  Fourier  path-integral  Monte 
Carlo.  A  genuine  Berezinskii-Kosterlitz-Thouless  transition  is  found  up  to  a  threshold  value 
g*  of  the  quantum  coupling,  beyond  which  no  phase  coherence  is  established.  Slightly  below 
g*  the  phase  stiffness  shows  a  reentrant  behavior  with  temperature,  in  connection  with  a  low- 
temperature  disappearance  of  the  superconducting  phase,  driven  by  strong  nonlinear  quantum 
fluctuations.  The  presence  of  dissipation,  due  to  the  shunt  resistences,  reduces  the  region 
where  this  reentrance  appears.  However  it  is  persisting  for  significant  dissipations  and 
eventually  vanishes  for  the  decreasing  of  the  shunt  resistences  when  the  quantum  fluctuations 
are  strongly  reduced  by  the  dissipation,  ects. 
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We  study  the  persistent  currents  of  normal/superconducting  loop  modulated  by  the  phase 
shift,  77,  of  electron  wave  functions  due  to  the  Andreev  reflection  process  at  the  boundaries. 
When  the  energy  splitting  E  is  smaller  than  the  pair  potential  A,  the  persistent  currents 
show  the  periodicity  of  the  superconducting  unit  flux  quantum,  O0  /  2,  regardless  of  the 

length  of  the  normal  and  the  superconducting  sector,  while  it  is  so  only  in  the  limiting  case 
of  thin  normal  segment  in  previous  studies.  In  the  limit  E  «  A ,  the  persistent  currents 

take  the  Josephson  type  formula,  1  <x  Va  sin(/r  /  2  -  2 rj) . 


Quantum  Order  in  Josephson  Junction  Arrays 

G.  Cristofano1,  V.  Marotta1,  A.  Naddeo2.  P.  Sodano3 4 

1  Dipartimento  di  Scienze  Fisiche,  Universita  di  Napoli  “Federico  II”  and  INFN, 
Sezione  di  Napoli,  Via  Cinthia  -  C.  U.  Monte  S.  Angelo  -  80126  Napoli ,  Italy 
2 Dipartimento  di  Scienze  Fisiche,  Universita  di  Napoli  “Federico  II”  and  Coherentia 
-  INFM,  Unita  di  Napoli,  Via  Cinthia  -  C.  U.  Monte  S.  Angelo  -  80126  Napoli,  Italy 
3 Dipartimento  di  Fisica,  Universita  di  Perugia  and  INFN,  Sezione  di  Perugia,  Via  A. 
Pascoli  -  06123  Perugia,  Italy 


Quantum  order  is  a  desirable  property  of  Josephson  junction  networks  since  it  may 
allow  for  their  use  as  devices  for  quantum  computation.  There  have  been  by  now 
several  proposals  [1],  In  this  contribution  we  show  that  folly  frustrated  Josephson 
networks  may  support  quantum  order  [2].  In  order  to  meet  such  a  request  we  analyze 
the  folly  frustrated  XY  (FFXY)  model  on  a  square  lattice  with  f/(l)®Z2  symmetry 
[3]  within  a  twisted  conformal  field  theory  (CFT)  framework  [4],  In  particular  we 
show  that  its  ground  state  is  degenerate,  the  different  states  being  accessible  by 
adiabatic  flux  change  techniques.  Such  a  degeneracy  is  shown  to  be  strictly  related  to 
the  presence  in  the  spectrum  of  quasiparticles  with  non  abelian  statistics  and  can  be 
lifted  non  perturbatively  through  vortices  tunneling,  so  giving  rise  to  a  highly 
entangled  state. 
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A  quantum  system  (S)  in  interaction  with  a  repeatedly  measured  one  (M) 
undergoes  a  non-unitary  time  evolution  pushing  it  into  some  specific  subspaces.  It  is 
assumed  that  the  measured  system  is  found  in  the  initial  state  at  each  step,  and  that  both 
the  M-S  interaction  and  its  duration  between  two  measurement  acts  are  the  same  during 
all  the  process.  Such  three  ingredients  (measured  state,  interaction,  interaction  duration) 
give  raise  to  a  “selection  rule”  determining  what  states  of  S  decay  (the  “residual” 
states)  and  what  states  are  decay-preserved  (the  “distilled”  states). 

On  the  basis  of  this  strategy,  we  report  a  method  for  distilling  the  vibrational  state 
of  a  trapped  ion  centre  of  mass  in  order  to  obtain  angular  momentum  eigenstates  and 
angular  momentum  Schrodinger  cat  superpositions. 
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The  quantum  dynamics  in  solid-state  nanostructures  has  attracted  much  attention  for 
building  blocks  of  quantum  computing.  We  present  some  recent  results  on  coherency, 
dissipation  and  dephasing  of  charge  and  spin  degrees  of  freedom  in  semiconductor  quantum 
dots. 

First,  we  discuss  coherency  of  charge  states  in  a  double  quantum  dot.  Consider  a  two-level 
system  consisting  of  two  charge  state,  in  which  an  excess  electron  occupies  an  energy  state  in 
one  or  the  other  dot.  When  these  two  states  are  coupled  by  a  tunneling  barrier,  one  can  induce 
coherent  charge  oscillations  by  means  of  non-adiabatic  excitation  with  a  high-speed  voltage 
pulse.  We  have  successfully  observed  clear  coherent  charge  oscillations.  The  charge  qubit  can 
be  manipulated  arbitrarily  by  tailoring  the  pulse  waveform.  The  experiments  indicate  the 
potential  application  for  quantum  information  processing. 

Second,  we  discuss  about  dissipation  of  spin  degree  of  freedom  in  a  quantum  dot.  Electron 
spin  in  quantum  dots  is  expected  to  have  a  long  decoherence  time.  We  investigate  the  energy 
relaxation  time  from  an  excited  state  to  the  ground  state  in  one-  and  two-electron  quantum 
dots.  We  find  that  the  relaxation  time  from  spin-triplet  state  to  spin-singlet  state  in  a 
two-electron  quantum  dot  can  be  longer  than  hundreds  of  microseconds.  The  long  spin 
relaxation  time  is  advantageous  for  spin  based  information  storage. 
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Extremely  weak  interaction  of  nuclear  spins  to  their  environments  yields,  on  one  hand,  their 
exceedingly  long  de-coherence  time,  but  on  the  other  hand,  makes  their  control  and  detection  a 
nontrivial  task.  Contriving  ingenious  means  of  controlling/detecting  nuclear  spins  is  thus  an  important 
experimental  challenge  for  exploring  the  implementation  of  spintronics,  quantum  information 
processing  or  quantum  computation  in  solid  state  devices.  We  show  here  that  spin-split  edge 
channels  in  the  integer  quantum  Hall  effect  (IQHE)  devices  provide  a  unique  tool  to  (i) 
polarize,  (ii)  unitary-transform,  and  (iii)  detect  nuclear  spins.  The  number  of  manipulated 
nuclear  spins  is  small,  on  the  order  of  109. 

The  one-electron  energy  of  a  finite  two-dimensional  electron  gas  (2DEG)  system  is 
quantized  into  discrete  Landau  levels  (LLs)  in  strong  magnetic  fields,  where  each  orbital 
LL  is  spin  split  into  spin-up  and  spin-down  Zeeman  sub-levels.  Due  to  confining 
potential,  the  LLs  increase  their  energies  as  a  boundary  of  the  2DEG  layer  is  approached 
and  form  respective  edge  channels  along  the  boundary  at  the  Fermi  level.  By  selectively 
transmitting  different  edge  channels  by  means  of  cross-gate  biasing  technique,  one  can 
unequally  populate  the  spin-up  edge  channel  and  the  spin-down  edge  channel  of  the 
lowest  orbital  LL  in  a  GaAs/AlGaAs  Hall  bar  device  in  the  v=2  IQHE  regime. 
Over-populated  edge  channel  (either  with  spin-up  or  spin-down)  can  be  chosen  by  the 
polarity  of  the  source-drain  voltage  applied  to  the  device.  When  the  up-spin  channel  is 
over  populated,  for  instance,  spin-up  electrons  undergo  spin-flip  transition  to  empty 
states  of  the  spin-down  edge  channel.  Due  to  the  hyperfine  interaction  [1]  this  process  is 
accompanied  by  the  down-to-up  spin-flop  transition  in  the  nuclear  spin  system,  leading  to 
the  development  of  strong  dynamic  nuclear  polarization  (DNP)  in  the  region  only  along 
the  edge  channels.  The  DNP,  in  turn,  induces  an  effective  magnetic  field  for  the  electron 
spin  system,  again,  via  the  hyperfine  interaction.  This  leads  to  a  change  of  the 
resistance  of  the  device,  through  which  the  DNP  is  sensitively  detected.  By  passing 
radio-frequency  (rf)  currents  through  a  fine  metal  strip  placed  near  the  edge  channels 
locally  generate  rf-magnetic  fields  in  the  region  of  DNP.  Nuclear  magnetic  resonance 
(NMR)  as  well  as  the  Rabi  oscillation  of  initialized  nuclear  spins  are  electrically  detected 


(through  the  resistance)  [2,3],  Free-induction  decay  as  well  as  the  spin-echo  experiments 
have  been  made,  suggesting  T2  =  80  ps.  In  addition,  by  externally  modifying  the 
confining  potential  via  side-gate-biasing  technique,  the  edge  channels  are  shifted 
side-wise  over  a  120nm  distance,  thereby  locally  probing  the  spatial  profile  of  the  DNP. 
DNP  proved  to  spread  roughly  over  a  20  nm-wide  region  (twice  the  magnetic  length). 
The  number  of  manipulated  nuclear  spins  is  suggested  to  be  on  the  order  of  109  The 
spin-split  edge  channels  are  thus  a  powerful  tool  to  (i)  efficiently  initialize,  (ii) 
unitary-transform,  and  (iii)  sensitively  detect  nuclear  spins. 
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The  Quantronium  [1]  is  a  superconducting  quantum  bit  circuit  (Fig.  1)  based  on  small 
Josephson  junctions  in  which  decoupling  from  the  external  circuit  ensures  a  long  coherence 
time.  The  state  of  the  qubit  is  manipulated  using  resonant  microwave  pulses,  or  adiabatic 
pulses  on  the  control  parameters  of  the  circuit  Hamiltonian.  We  demonstrate  here  that  NMR 
methods  can  be  efficiently  used  for  manipulating  the  qubit  and  probing  decoherence 
mechanisms  [2], 

First,  we  have  shown  that  arbitrary  transformations  can  be  implemented  by  combining 
microwave  pulses.  Two  pulse  sequences  with  different  phases,  that  implement  rotations 
around  orthogonal  axes,  have  been  performed.  The  overall  agreement  with  theory  shows  that 
rotations  combine  as  expected  (Fig;  2).  We  have  addressed  the  issue  of  transformation 
robustness  using  the  composite  pulses  developed  in  NMR.  We  have  shown  in  particular  that 
the  CORPSE  sequence  efficiently  compensates  for  of  resonance  effects. 

In  order  to  investigate  decoherence  mechanisms  acting  in  the  quantronium,  we  have 
developed  new  methods  to  probe  quantum  coherence.  In  particular,  we  have  used  two  pulse 
experiments  with  a  temporary  change  of  the  working  point  during  the  time  delay  between  the 
two  pulses,  and  implemented  echo  experiments  (see  Fig  3). 

We  have  also  measured  decoherence  during  driven  evolution  in  presence  of  microwaves. 
We  have  measured  the  decay  of  Rabi  oscillations,  and  of  the  spin-locking  polarisation 
following  the  usual  spin-locking  protocol  in  NMR  (see  Fig  3).  We  find  that  the  decay  time  of 
Rabi  oscillations  and  of  the  spin-locking  polarisation  can  be  longer  than  during  free  evolution 
coherence  time. 

These  driven  regime  and  spin  echo  experiments  rise  the  question  of  increasing  the  effective 
coherence  time;  which  is  a  central  issue  for  qubits. 


Fig  1:  top:  circuit  diagram  of  the  quantronium.  The 
Hamiltonian  is  controlled  by  two  parameters:  the 
charge  coupled  to  the  box  island  between  the  two 
small  junctions,  and  the  phase  across  their  series 
combination.  The  qubit  state  is  manipulated  by 
applying  resonant  microwave  pulses,  or  adiabatic 
pulses  on  the  control  parameters.  The  oscillations  of 
the  switching  probability  of  the  readout  junction  reveal 
the  Rabi  oscillations  of  the  qubit  state. 


[1]  D.  Vion  et  al..  Science  296  (2002). 
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Fig.2:  Switching  probability  after  two  %!2  pulses 
with  different  phases,  corresponding  to  different 
rotation  axes.  The  Ramsey  patterns  are  shifted  as 
predicted  for  each  phase  combination.  This 
experiment  proves  that  rotations  can  be 
combined  in  order  to  implement  any  qubit 
transformation. 
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Fig.  3:  Top  :  switching  probability  (dots)  after  a 
Ramsey  {tc/2(X),  it/2(X)}  sequence  as  a  function 
of  the  time  delay  between  pulses.  Fits  of  the 
envelope  (  time  constant  350  ns).  Middle: 
example  of  echo  measured  in  a  {jt/2(X),  n(X), 
jc/2(X)}  sequence.  Thin  line:  echo  signal  at  the 
nominal  minimum  position,  bold  line:  exponential 
fit  of  the  envelope  (  550  ns  time  constant), 
dashed  line:  fit  of  the  lower  envelope  of  the 
Ramsey  pattern  measured  in  the  same  conditions 
(220  ns  time  constant).  Bottom:  switching 
probability  after  two  spin-locking  sequences,  as  a 
function  of  the  sequence  duration.  Thick  lines: 
exponential  fits  of  the  envelopes,  with  time 
constant  650  ns.  The  dashed  lines  show  a  fit  of 
the  envelope  of  the  Ramsey  pattern  measured  in 
the  same  conditions  (time  constant :  320  ns). 
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The  Quantronium  [1]  is  a  superconducting  Josephson  quantum  bit  (Fig.  1)  whose  quantum 
state  can  be  arbitrarily  transformed  by  combining  quasi-resonant  microwave  pulses  and 
adiabatic  pulses  [2]  on  two  control  parameters  of  the  circuit  Hamiltonian:  an  electric  charge 
Ng  and  a  phase  8. 

The  coherence  time  T2  of  this  qubit  during  its  free  evolution  was  measured  by  several 
different  pulse  sequences,  and  at  different  working  points  (Fig.  2)  where  the  dominant 


contribution  to  decoherence  is  the  external  charge  or  phase  noise  at  first  or  second  order.  The 
coherence  time  during  driven  evolution  was  also  measured  as  a  function  of  the  detuning 
between  the  microwave  frequency  and  the  qubit  frequency  (Fig.  3).  All  these  data  are 
analyzed  using  a  generalized  Bloch-Redfield  model. 


[1]  D.  Vion  et  al.,  Science  296  (2002). 

[2]  E.  Collin  et  al.,  cond-mat/0404503. 


Fig  1:  Circuit  diagram  of  a  qubit 
prototype:  the  quantronium.  The  qubit 
state  is  manipulated  by  applying 
quasi-resonant  microwave  pulses  on 
Ng=CgU/(2e),  and  adiabatic  pulses  on 
Ng  or  Ib.  The  state  is  read  out  by  pulsing 
Ib  to  the  critical  current  of  the  large 
Josephson  junction. 
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Fig.2:  Coherence  times  T2  measured  by  different 
techniques  as  a  function  of  the  phase  or  charge 
shift  from  the  optimal  working  point  (8  =0  ; 
Ng=T/2).  Lines  correspond  to  the  model. 
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Fig.3:  Coherence  times  during  microwave 
driven  Rabi  oscillations  as  a  function  of 
the  detuning.  Line  corresponds  to  the 
model. 
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Although  Josephson  junction  quantum  bits  show  great  promise  for  quantum  computing, 
the  origin  of  dominant  decoherence  mechanisms  remains  unknown.  All  superconducting 
qubit  experiments  thus  far  have  shown  relatively  low  (coherence)  amplitude  and  short 
relaxation  times  (for  Rabi  oscillations).  Recently,  we  have  improved  the  operation  of  a 
Josephson  Junction  based  "phase"  qubit  by  incorporating  it  into  a  superconducting  loop. 
This  has  allowed  us  to  further  decouple  the  qubit  from  its  environment  and  has 
eliminated  quasiparticle  heating  effects.  A  new  detection  scheme  based  on  an 
asymmetric  d.c.  SQUID  allows  detection  of  the  quantum  states  with  high  efficiency  and 
a  controllable  back  action.  This  new  device  has  revealed  previously  unknown  two  state 
microwave  resonators  within  the  qubit  Josephson  junction  itself. [1]  It  seems  likely  that 
all  Josephson  junction  based  qubits  are  suffering  from  these  material  defects  as  well. 
Although  two-level  systems  are  known  to  exist  in  amorphous  materials[2],  these 
individual  two-level  microwave  resonators  within  Josephson  tunnel  barriers  have 
remained  hidden  for  over  forty  years.  As  major  sources  of  decoherence,  they  play  a  vital 
role  in  understanding  the  full  effective  Hamiltonian  of  Josephson  tunnel  junctions. 

At  certain  bias  points,  strong  coupling  to  these  resonators  destroys  the  coherence  of  the 
qubit.  Although  previous  reports  have  focused  on  the  characteristic  decay  time  of 
coherent  oscillations,  our  data  demonstrates  that  decoherence  from  these  spurious 
resonators  primarily  affects  the  amplitude  of  the  oscillations.  We  have  constructed  a 
model  for  the  microscopic  resonators,  where  we  consider  two-level  states  in  the  barrier 
with  large  tunneling  matrix  elements  corresponding  to  a  microwave  frequency.  These 
two  level  systems  couple  to  the  qubit  by  producing  critical  current  fluctuations  in  the 
tunnel  barrier.  If  we  assume  that  the  two  states  of  the  resonator  correspond  to  ON  and 
OFF  configurations  of  individual  conduction  channels  in  the  tunnel  barrier,  then  their 
strength  can  be  estimated  from  the  current-voltage  characteristics  of  the  junction. 
Essentially,  low  subgap  currents  indicate  weaker  resonators  and  improved  qubit 
performance.  This  suggests  that  the  quality  of  the  Josephson  junction  determines  the 
density  and  size  of  these  resonators  and  therefore  the  quality  of  the  qubit. 

State-of-the-art  Josephson  junctions  employing  superconducting  (A1  or  Nb)  electrodes 
and  native  oxide  tunnel  barriers  typically  are  fabricated  using  sputter  deposition  onto 
thermally  oxidized  Si  wafers  for  microelectronics  fabrication.  For  most  applications, 
such  as  SQUID  technologies,  there  has  never  been  a  necessity  to  further  improve  the 
crystalline  quality  of  the  junctions  to  correspondingly  improve  the  device  performance. 
We  believe  that  these  spurious  resonators  originate  from  bulk  or  interface  defects  and 
inhomogeneous  or  non-uniform  oxidation  during  tunnel  barrier  formation.  In  this  work, 
we  detail  recent  efforts  to  engineer  Josephson  junctions  with  increasingly  improved 


crystalline  quality  and  smoother  interfaces,  and  to  correlate  the  microstructure  of  the 
junctions  with  low  frequency  transport  measurements  and  quantum  measurements  of 
phase  qubits.  Our  data,  thus  far,  indicates  a  strong  correlation  between  improved 
“smoothness”  of  the  superconductor-insulator  interface  and  junction  performance.  By 
creating  tunnel  junctions  from  sputtered  trilayers,  as  opposed  to  a  two-step  “ion-milled” 
process,  we  have  reduced  the  size  of  subgap  currents  and  have  shown  improved 
coherence  of  fabricated  phase  qubits. 

We  will  discuss  the  details  of  samples  grown  using  two  techniques,  sputtered  trilayers 
and  evaporated  trilayers.  Sputtered  trilayers  are  a  simple  to  fabricate  and  are  currently 
the  preferred  method  for  qubit  fabrication.  The  evaporated  trilayers,  on  the  other  hand, 
are  grown  using  a  truly  unique  method,  and  while  in  its  preliminary  stages,  this  growth 
technique  may  someday  replace  sputtered  trilayers  for  qubit  applications.  In  order  to 
prepare  the  evaporated  trilayers,  we  have  grown  tunnel  junctions  using  an  atomically 
flat  metallic  seed  layer  of  A1  on  Si(lll)-(7  x  7)  substrates.  The  seed  layers  were 
deposited  at  120  K  and  subsequently  annealed  at  425  K  to  achieve  ultra  smooth, 
epitaxial  films.  The  Josephson  trilayer  was  completed  through  the  homoepitaxy  of  A1 
deposited  on  the  seed  layer  at  room  temperature,  oxidized  in  a  conventional  manner  and 
then  capped  with  the  A1  counterelectrode  in  situ  under  ultra  high  vacuum  conditions. 
These  evaporated  trilayers  were  fabricated  into  tunnel  junction  devices  of  various 
diameters  and  transport  measurements  were  made  at  ~  30  mK.  I-V  characteristics  for 
samples  with  various  degrees  of  crystalline  quality  have  been  compared,  and  we 
conclude  that  there  is  a  strong  correlation  between  the  morphology  of  oxidized  base 
electrodes  that  are  smooth  and  the  lowering  of  subgap  currents.  This  is  an  important 
step  in  addressing  the  origin  of  spurious  resonators  ubiquitously  found  in  our  Josephson 
phase  qubits. 

In  order  to  eliminate  the  effect  of  these  resonances  completely,  further  materials 
research  is  being  investigated.  We  have  made  more  epitaxial  trilayers,  added  impurities 
to  the  superconducting  leads,  as  well  as  used  various  oxidation  methods.  We  will 
present  our  preliminary  results.  Furthermore,  simple  changes  have  been  made  in  order 
to  investigate  the  ways  in  which  these  resonators  can  influence  junction  performance. 
We  have  changed  the  area  of  the  qubit  junction,  the  inductance  of  the  superconducting 
loop,  and  the  mutual  inductance  to  its  bias  coil.  We  have  confirmed  that  a  reduction  in 
area  does  reduce  the  number  of  resonators  within  the  frequency  range  of  operation,  but 
each  resonator  has  an  increased  strength.  We  have  measured  the  distribution  of  these 
resonators  in  terms  of  their  strength  and  their  number  density  over  the  measurable 
frequency  range.  This  data  helps  to  further  characterize  junction  quality.  We  have  also 
developed  a  mathematical  model  of  these  resonators  which  agrees  well  with 
experimental  data.  This  model  provides  insight  both  into  the  possible  physical  origins  of 
the  resonances,  as  well  as  into  refocusing  schemes  to  eliminate  their  effects  on  the  qubit. 

[1]  R.  W.  Simmonds,  K.  M.  Lang,  D.  A.  Hite,  D.  P.  Pappas,  and  John  M.  Martinis,  cond-mat/0402470 
(2004) 
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Here  we  describe  a  novel  motor,  which  is  capable  of  non-hysteretic  (reversible) 
positioning  of  a  sample  in  three  orthogonal  directions  (x,  y,  and  z)  with  nanometer 
resolution  both  at  room  temperature  and  at  cryogenic  temperatures.  The  motor  can 
translate  the  sample  ±  1.5  mm  in  the  x-y  plane  and  ±  2mm  in  the  z-direction  and 
within  this  range  park  the  sample  at  any  position  in  firm  contact  with  the  motor  body 
(zero  voltage  on  the  piezo  stacks). 

The  motor  is  integrated  with  a  small  specially  designed  head,  which  except  for 
modifications  necessary  for  cryogenic  operation  assembles  a  standard  insert  for 
Scanning  Probe  Microscopy  (SPM).  The  SPM  head  faces  the  sample  stage  and 
consists  of  a  conventional  piezoelectric  x-y-z  tube  scanner  carrying  either  an  STM  tip 
(STM-head)  or  an  AFM  cantilever  (AFM-head).  With  the  sample  parked  in  fixed 
position  by  the  x-y-z  motor  the  SPM  head  allows  for  fine  positioning  and  scanning 
with  atomic  precision.  A  new  millikelvin  AFM  head  is  under  construction.  It  utilizes 
single  mode  glass  fibers  both  for  supplying  the  laser  beam  focused  onto  the  cantilever, 
and  for  transmitting  the  reflected  beam  back  to  the  room  temperature  detection 
system. 

The  x-y  translation  of  the  sample  is  facilitated  by  two  sets  (A  and  B)  of  piezo  stacks 
(“legs”)  glued  onto  a  thick  brass  plate.  The  legs  are  placed  evenly  (60  degrees)  in  a 
circle  with  diameter  approx.  20mm  (see  Fig.  1).  Each  leg  consists  of  three  5mm  by 
5mm  piezo  stacks  (“legs”),  each  with  four  x-shear  plates  glued  in  series  with  four  y- 
shear  plates  and  a  multilayer  z-stack.  On  top  and  bottom  of  each  leg  is  glued  a  5mm 
by  5mm  A1203  plate.  After  gluing  the  six  legs  are  optically  polished  to  obtain  same 
length.  The  operational  relies  on  the  well  known  louse  principle,  modifyed  to  a  6  leg 
configuration,  where  in  every  step  3  legs  are  lifted  fully  free  from  the  polished  quartz 
or  sapphire  surface  before  they  are  jointly  moved  to  the  new  position.  The  z 
translation  utilizes  a  similar  working  principle  that  enables  a  linear  motion. 

The  motor  can  be  used  as  a  coarse  positioning  device,  which  with  a  (large)  number  of 
(up  to  a)  micrometer  long  x-y-z  steps  rapidly  brings  the  selected  region  of  the  sample 
in  first  contact  with  the  probe  tip  in  the  SPM  head.  Contact  is  registered  as  either  the 
first  small  tunnel  current  in  the  STM  head  or  the  small  change  in  the  deflection  of  the 
laser  beam  focused  on  the  cantilever  in  the  AFM  head.  Maintaining  the  contact  the 
motor  then  parks  the  sample  with  zero  voltage  on  the  motor  piezos.  Here  the  galvanic 
contacts  to  the  wires  leading  to  the  motor  piezos  may  be  disconnected,  so  that  the 
electrical  noise  from  the  high-tension  amplifiers  cannot  induce  mechanical 
disturbances  in  the  sample  position.  After  having  parked  the  sample  x-y-z  scanner 
tube  unit  regulates  the  position  of  the  STM  tip  or  AFM  cantilever.  Actually,  one  may 
use  very  small  x-y-z  steps  in  the  motor  to  adjust  the  tip-sample  distance  so  that  also 


the  off-set  voltages  applied  to  the  tube  scanner  are  near  zero.  Recent  experimental 
results  will  be  reported  on. 

Presently  we  are  investigating  the  creep  of  the  motor  piezos  and  the  tube  scanner  at 
low  temperatures.  For  the  non-scanning  applications  eg.  nano-lithography  and  nano¬ 
junction  fabrication  it  is  necessary  to  have  a  linear  dependence  of  the  tip  position  on 
the  piezo  voltage.  However,  noise  in  the  “feed  back”  optical  loop  used  is  estimated  to 
induce  a  1-10  nanometer  fluctuation  in  the  spatial  resolution.  In  many  cases  this  limits 
the  accuracy  and  for  better  precision  one  has  to  “live  with”  the  creep  in  the  piezos. 
The  optical  servo  system  will  be  briefly  described. 

Together  the  motor  and  the  SPM  head  form  an  instrument,  which  enables  atomic 
resolution  both  in  the  scanning  modes  and  at  fixed  positions.  This  allows  for  a  wide 
spectrum  of  applications,  such  as  nano-lithography  and  studies  involving  in  situ 
adjustable  junctions.  A  particular  advantage  of  cryogenic  operation  is  that  thermal 
fluctuations  in  the  piezo-electrical  materials  are  very  small.  The  disadvantage  is  that 
the  piezo-electrical  constants  at  e.g.  4.2K  are  reduced  to  about  20%  of  the  room 
temperature  values.  The  motor  mounted  with  an  SPM  head  is  very  compact  and  fits 
into  the  34mm  diameter  bore  of  a  standard  (Oxford)  millikelvin  3He/4He  cooler 
placed  in  the  center  of  a  10  tesla  superconducting  magnet. 

The  motor  and  the  SPM  head  are  jointly  run  by  a  computer  (with  an  A/D  interface 
input/output  card)  in  conjunction  with  a  standard  SPM  controller  (Veeco,  Digital 
Instruments).  Depending  on  the  controller  and  the  software  one  may  utilize  all  the 
standard  modes  of  operation,  eg.  tapping  mode. 


